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ABSTRACT 


A  code  system  (called  UNC-SAM)  which  evaluates  the  transport  of  neutrons  or 
gamma  rayo  using  Monte  Carlo  methods  in  complex  three-dimensional  geometry, 
has  been  written  in  FORTRAN  for  the  Control  Data  Corporation  1604 -A  computer. 
The  code  will  calculate  fluxes,  flux  dependent  functionals  such  us  doses,  and  their 
standard  deviations  in  geometry  comprised  of  rectangular  parallelepipeds  which, 
in  turn,  may  contain  spheres,  cylinders,  parallelepipeds,  or  wedges. 

In  order  to  increase  the  efficiency  of  the  Monte  Carlo  transport  game,  UNC-SAM 
uses  an  importance  sampling  technique.  In  addition,  in  the  particular  case  of 
evaluating  neutron  fluxes  in  small  volume  detectors,  a  scoring  by  analytical  esti¬ 
mation  referred  to  as  “flux  at  a  point”  is  employed 
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1.  INTRODUCTION 


The  United  Nu:;lear  Corporation  Stochastic  Approximation  Method  (UNC-SAM)  is 
a  code  system  in  the  FORTRAN  language  for  evaluating  the  transport  of  neutrons 
or  photons  through  matter  by  Monte  Carlo  techniques.  It  handles  a  three-dimen¬ 
sional  geometry  composed  of  rectangular  parallelepipeds  which,  in  turn,  may  con¬ 
tain  spheres,  cylinders,  parallelepipeds,  or  wedges. 

UNC-SAM  represents  the  latest  result  of  a  continuing  effort  to  develop  complex 
geometry  Monte  Carlo  programs  of  ever  greater  versatility  and  speed.  *  Most 
practical  reactor  and  shield  geometries  can  be  handled  by  the  geometry  routines. 
Despite  this,  the  geometry  tracking  portions  of  a  problem  consume  considerably 
less  than  half  the  problem  running  time.  In  addition  to  high  efficiency  geometry 
routines,  region  and  energy  weighting  and  a  special  “flux-at-a-point”  routine 
have  been  incorporated  to  reduce  problem  running  time  when  fluxes  are  required 
at  improbable  regions  of  phase  space. 

Some  applications  in  which  UNC-SAM  (or  its  predecessor  ADONIS^)  have  been 
used  successfully  are: 

1.  Transmission  of  neutrons  through  ducted  shields 

2.  The  analysis  of  experiments  aimed  at  determining  the  production  of 
gamma  rays  in  shield  materials  by  high  energy  neutrons 

*Most  of  the  numerical  techniques  employed  in  UNC-SAM  were  either  formulated 
or  invented  by  M.  H.  Kalos  who  guided  the  program  throughout  its  development. 
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3.  The  transmission  of  neutrons  and  gamma  rays  through  air,  taking  account 
of  air  and  ground  scattering 

4.  Investigation  of  shielding  provided  by  proposed  reactor  containment  struc¬ 
ture 

5.  Measurement  of  shielding  effocuveness  of  nested  boxes  containing  geo¬ 
metrically  complex  devices 

6.  Weapon  vulnerability  studies. 

The  code  system  is  written  in  FORTRAN  for  the  Control  Data  Corporation’s  (CDC) 
1604 -A  computer.  A  32  K  core  and  <  6  magnetic  tapes  are  required.  UNC-SAM 
is  in  reality  a  chained  series  of  independent  programs  which  process  cross-sec¬ 
tion  and  geometry  data,  do  the  transport  problem,  and  edit  the  results.  A  brief 
summary  of  the  independent  programs  contained  in  UNC-SAM  follows.  Fig.  1  in¬ 
dicates  the  information  flow  in  the  system. 

GENDA 

The  Generation  of  Data  routine  (GENDA)  starts  with  a  fundamental  set  of  neutron 
or  photon  cross-section  data  tabulated  as  a  function  of  energy.  It  allows  for  the 
basic  tabulation  of  the  total  cross  section,  the  scattering  cross  secfon,  the 
Legendre  expansion  coefficients  of  the  differential  scattering  cross  section,  the 
inelastic  scattering  cross  section  for  both  level  and  continuum  scattering,  etc. 
These  data  are  tabulated  generally  on  differing  energy  meshes.  GENDA  processes 
the  data  to  generate  a  set  of  cross-section  data  on  a  desired  energy  mesh.  It  in¬ 
cludes  the  option  of  averaging  with  respect  to  a  given  weighting  function,  pre¬ 
scribed  by  input,  if  desired.  Otherwise,  it  determines  the  desired  cross-section 
set  by  either  linear  interpolation  or  by  log-log  interpolation.  The  spectrum  of 
neutrons  emitted  by  inelastic  scattering,  according  to  statistical  theory,  is  used 
for  heavier  materials.  For  lighter  materials,  such  as  Li  and  Be,  special  treat¬ 
ments  are  available.  The  output  of  this  routine  is  a  tape  which  may  be  used  as 
input  to  GENPRO. 
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GENPRO 


The  GENDA  Processor  (GENPRO)  routine  uses  the  GENOA  output  to  determine 
the  tables  of  probabilities  peculiar  to  a  nuclide,  required  for  a  Monte  Carlo  code, 
at  a  prescribed  energy  mesh.  For  instance,  the  probability  that  a  neutron  suffers 
an  elastic  scattering  is  given  by  the  ratio  of  elastic  to  total  scattering  at  that  en¬ 
ergy.  The  probability  of  inelastic  scattering  from  one  energy  to  another  is  com¬ 
puted  in  this  routine.  At  its  conclusion,  one  has  a  library  tape  called  the  Element 
Data  Tape  (EDT),  on  a  given  energy  mesh  to  be  used  in  subsequent  processing. 

Up  to  this  point,  no  problem  oriented  input  was  needed.  The  functioning  of  later 
routines  (DATORG,  EZGEOM,  VANGEN,  GASP)  will  be  dependent  on  such  problem 
oriented  data. 

DATORG 

Given  the  concentrations  of  the  materials  comprising  the  configuration,  an  orga¬ 
nized  data  tape,  bearing  the  total  macroscopic  cross  sections  as  well  as  other 
probabilities  for  mixtures  of  nuclides,  is  created  by  the  DATORG  routine.  These 
are  the  problem -dependent  cross-section  data  tabulated  against  energy  at  the 
prescribed  energy  mesh  previously  used  in  GENDA  and  GENPRO. 

EZGEOM 

This  is  a  routine  which  takes  a  simplified  geometrical  description  of  the  physical 
system,  as  provided  by  the  problem  originator,  and  produces  the  rather  complex 
set  of  data  required  by  the  transport  program  ADONIS. 

The  problem  originator  must  decompose  space  into  a  set  of  boxes  known  as  “or¬ 
dinary  regions.”  Inside  these  boxes  may  be  placed  “nonordinary  regions”  con¬ 
sisting  of  spheres,  cylinders,  wedges,  or  other  boxes.  Nonordinary  regions  may 
enclose  other  nonordinary  regions. 
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VANGEN 


Source  input  corresponding  to  a  primary  neutron  or  gamma -ray  problem  is  pro¬ 
vided  through  the  Volume  Anisotropic  Source  Generator  (VANGEN).  VANGEN 
creates  a  source  tape  with  the  following  types  of  distributions: 

P-(i)  point  source 

P-(ii)  uniform  plane  source 

P-(iii)  uniform  volume  source. 

Corresponding  to  each  of  the  above,  a  choice  of  angular  distributions  describing 
the  source  is  available. 

A-(i)  isotropic 

A-(ii)  anisotropic  according  to  a  given  (input)  distribution 

A-(iii)  monodirectional 

A-(iv)  isotropic  in  the  half  plane. 

Finally,  to  describe  the  energy  variation  of  the  source,  these  choices  are  available: 

E-(i)  fission  spectrum 
E-(ii)  monoenergetic 

E-(iii)  an  arbitrary  distributed  energy  source  with  equal  probability  steps 
prescribed  by  input. 

The  output  of  VANGEN  is  a  tape  for  further  processing. 

GASP 

The  routine  Gamma  Source  Particle  generator  (GASP)  computes  gamma -ray  pro¬ 
duction  from  a  primary  neutron  problem  to  serve  as  a  source  for  a  secondary 
gamma-ray  problem.  GASP  requires: 

1.  An  interaction  tape  (generated  by  a  previous  ADONIS  problem)  containing 
a  record  of  all  interactions  able  to  cause  gamma-ray  emission 
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2.  Probability  numbers  for  each  nuclide  of  the  problem  giving  the  probable 
number  of  gamma  rays  emitted  as  a  function  of  energy 

3.  An  angular  distribution  (as  in  VANGEN)  for  newly  created  gamma  rays. 

The  routine  can  cope  with  the  generation  of  other  secondaries  upon  modification 
of  the  interaction  tape  input.  The  output  of  either  GASP  or  VANGEN  is  a  source 
particle  tape  which  becomes  input  to  ADONIS. 

ADONIS 

Given  the  appropriate  output  tapes  of  DATORG,  EZGEOM,  VANGEN,  or  GASP, 
and  parametric  input  defining  the  problem  to  be  solved,  one  is  now  ready  to  pro¬ 
ceed  with  the  neutron  or  gamma  transport  code  (ADONIS)  The  program  starts 
sets  of  100  initial  source  particles  from  the  source  tape  and  processes  these  par¬ 
ticles  through  all  of  their  interactions  and  migrations,  until  death,  reflection, 
transmission,  or  degradation  below  an  energy  cutoff  occurs.  As  the  particle  trav¬ 
els,  a  variety  of  information  is  recorded,  including: 

1.  An  interaction  tape  for  future  secondary  emission  processing,  or  appro¬ 
priate  response  function  calculations  relating  to  dose  calculations,  etc. 

2.  The  total  flux  as  a  function  of  energy  in  each  of  the  ordinary  and  non¬ 
ordinary  geometric  figures,  along  with  Its  standard  deviation 

3.  A  transmitted  particle  tape  for  future  processing. 

Some  of  the  Monte  Carlo  techniques  employed  in  the  ADONIS  code  are  listed  below. 

1.  Splitting  or  Russian  roulette  is  performed  at  each  geometric  boundary 
crossed  during  the  process  of  reaching  a  collision  point.  The  splitting 
process  is  determined  by  the  weights  assigned  to  the  regions  on  either 
side  of  the  boundary.  The  weight  of  any  region  may  vary  with  the  neutron 
energy-.  Any  additional  particles  generated  by  the  splitting  process  are 
referred  to  as  latent  particles. 
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2.  Fluxes  in  small  detector  regions  are  computed  by  the  method  referred  to 
as  flux-at-a-point.  (At  present  this  applies  to  neutron  problems.) 

Some  of  the  output  provided  by  ADONIS  includes: 

1.  Number  of  absorptions,  births,  deaths,  degradations 

2.  Flux  and  standard  deviation  per  region  energy  bin 

3.  Dose  and  standard  deviation  per  region. 


SAM 

SAM  is  the  main  program  of  an  overlay  tape  containing  the  DATORG,  EZGEOM, 
VANGEN,  GASP  and  ADONIS  programs.  SAM  functions  as  a  program  sequencer 
reading  the  title  card  of  each  set  of  input  for  a  given  problem  calling  the  program 
required.  It  allows  a  user  to  place  his  input  deck  in  any  manner  providing  that, 
when  a  program  appears  in  the  sequence,  ail  programs  that  provide  input  to  that 
program  have  already  been  run. 

The  following  sections  of  this  report  present  a  detailed  discussion  of  the  individu¬ 
al  routines  above.  Where  applicable,  the  physical  principles  underlying  the  pro¬ 
gram  are  discussed. 
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2.  THE  GENERATION  OF  DATA  ROUTINE  (GENDA) 

Given  experimental  or  theoretical  microscopic  cross-section  data  at  arbitrary 
energies  in  a  given  range  for  elements  of  interest,  it  is  often  desirable  to  process 
these  data  into  cross-section  sets  at  specified  energy  point  values. 

The  GENDA  routine  has  been  programmed  in  FORTRAN  to  accomplish  this  in  the 
following  manner: 

1.  The  raw  cross-section  data  are  tabulated  according  to  an  established 
format  which  can  be  easily  modified  and  updated.  The  point  values  of 
the  input  cross-section  data  are  selected  to  give  the  best  detailed  ’-epre- 
sentation  of  the  cross-section  variation  with  energy  for  the  interpolation 
technique  used. 

2.  For  a  given  element  and  type  of  cross  section,  interm.ediate  point  values 
are  obtained  from  these  tables,  using  either  linear  or  log-log  interpola¬ 
tion  techniques  as  directed  by  'he  input. 

3.  Weighted  point  values  of  cross  sections  in  a  specified  energy  mesh  may 
be  evaluated  using  a  given  weighting  function.  This  input  function  is 
tabulated  in  a  similar  manner  as  the  raw  cross-section  data  and  inter¬ 
mediate  point  values  are  obtained  by  interpolation. 


9 


2.1  TECHNICAL  DESCRIPTION 


Let  a(E)  denote  generically  a  t5T)ical  cross  section  at  energy  E.  This  might  stand 
for  the  total  cross  section,  the  scattering  cross  section,  the  Legendre  expansion 
coefficient  (fp)  of  the  differential  scattering  cross  section,  the  inelastic  scattering 
cross  section  for  level  or  continuum  scattering,  etc.  (the  cross-section  types  in¬ 
cluded  appear  in  the  input  instructions).  Let  aQ(E)  be  the  continuous  portion  of  the 
typical  cross  section  a{E).  On  occasion,  where  many  resonances  exist,  a  set  of 
resonance  parameters  Aj^(Ej^)  may  be  given  for  each  generic  cross  section  type 
such  that  for  the  purpose  of  suitable  averages,  a(E)  is  representable  as: 

•^max 

a(E)  =  ao(E)  +  ^  A^(E^)  6{E-E^)  (1) 

where  6(E-Ei/)  is  the  Dirac  delta  function  of  argument  E-Ej^.  This  characteriza¬ 
tion  of  cr(E)  is  allowed  only  when  the  cross-section  averaging  routine  (CAR)  is 
used;  otherwise,  an  infinitely  narrow  resonance  peak  is  not  physically  reasonable, 
nor  computationally  manageable. 

Let  Ej  ,  j  =  1,  2,  ...,  J  be  a  set  of  final  energies  in  decreasing  order  at  which  the 
output  data  are  to  be  tabulated.  These  may  be  tabulated  either  as  input  or  com¬ 
puted  given  the  highest  energy,  E^^ax*  the  number  of  final  energies,  J,  and  the 
constant  lethargy  spacing,  h,  such  that 

h  =  llnl^  (2) 

J  Ej 

Where  Ej  is  the  lowest  energy. 

Let  0(E)  be  an  arbitrary  weighting  function  supplied  by  an  input  table 
’’’max 

0(E)  =  0o(E)  +  J)  Br(E^)  6{E-Er).  (3) 

T=1 
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(The  E-j-  must  not  have  an  energy  value  equal  to  any  of  the  Ej^.)  Let  ag(E)  and  <Pq{E) 
be  given  as  functions  of  eneigy  such  that  between  tabulated  values  these  functions 
are  either  linear  or  exponential. 


The  GENDA  code  then  provides  for  the  following  mutually  exclusive  options. 


1.  Find  a(Ep);  j  =  1,  ...,  by  linear  interpolation 

2.  Find  cr(Ej  );  j  =1,  ...,  Jj^ax  log-log  interpolation 

3.  Find  a(Ep);  j  =  2,  ...,  J-1  by 


X 


(Ep+Ef.,)/2 


a(EG)  = 


0q(Ej^)  +  2Bj(E-^)  cTq(E7-) 


X 


(Ep+Ep_j)/2 
(Ep+  Ep+,)/2 


0o(E)  dE  +  2Br(E7.) 


(4) 


where  the  sums  include  those  values  of  Ej^  and  Ej  falling  within  the  integration 
interval. 


For  the  purpose  of  CAR,  the  energy  mesh  used  is  made  up  of  the  set  of  all  ener¬ 
gies  at  which  a  and  are  given,  plus  the  end  points  of  the  interval  of  integration. 
In  the  present  program,  the  integration  is  carried  out  numerically,  using  the 
trapezoidal  rule.  Consequently,  the  richer  the  input  mesh,  the  better  the  results. 
It  should  be  noted  that  exact  integration  could  be  performed  at  the  price  of  some 
modifications  of  the  program. 

The  cross  section  at  thermal  energy  may  be  included  in  GENDA.  In  that  case,  the 
lowest  entry  in  the  final  energy  mesh  tabulation  is  the  thermal  energy. 

In  addition,  GENDA  can  compute  the  synthesized  differential  scattering  cross  sec¬ 
tion 
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0(eP,u.)  = 


<^n(Ep) 

Av 


L 


S 

je=o 


(2X  +  1)  f£(Ef)  P£(aj) 


for  w  =  -1  (0.01)  1  (5) 


where  =  the  elastic  scattering  cross  section 

PjfCco)  =  the  Legendre  polynomial  of  order  £  evaluated  at  the  angle  cosine  cu 

p 

f^(Ej  )=  the  Legendre  expansion  coefficient. 

All  negative  values  for  a(Ep,a))  are  printed  out.  An  edit  on  the  coarser  mesh 
-1  (0. 1 )  1  is  also  available. 

For  the  treatment  of  inelastic  scattering,  it  is  necessary  to  describe  in  greater 

p 

detail  the  kernel  g(Ej  ,E')  which  represents  the  probability  that,  following  an  in¬ 
elastic  scattering,  a  neutron  at  energy  Ep  will  scatter  into  energy  E'.  The  spacing 
of  the  discrete  energy  levels  in  the  nucleus  is  known  to  decrease  with  increasing 
energy.  At  high  energies,  the  levels  are  so  close  that  they  may  be  considered 
as  a  continuum  of  energy  levels.  Let  EgQ  be  the  energy  at  which  the  discrete 
level  spectrum  can  be  considered  to  pass  over  into  a  continuous  spectrum.  Then, 
for  incident  neutron  energies  less  than  or  equal  to  Eg(j,  only  discrete  level  in¬ 
elastic  scattering  occurs.  The  kernel  g(Ep,E')  is  representable  as  a  sum  of 
delta  functions 


g{E^,E')  =  gjy{Ef,E'}  =  2  MEP)  ^[Ef-CE'  +  E^^)]  E^<E^^E^^ 

all  u 

(6) 

where  =  the  energy  level 

E^  =  the  energy  of  the  discrete  level 
aj^(Ep)  =  the  probability  that  a  neutron  at  energy  Ep  will  scatter  into  energy 
E'  =  Ep  -  E^ 
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Fi'om  the  definition  of  gj)  it  follows  that 


gD(Ef,E')  dE'  = 


E  a^{Ep)  =  1 

all  1/ 


(7) 


for  Ep  5^  E'.  If  Ep  <  E',  the  inelastic  scattering  cross  section  is  zero  and 
gD(Ep,E0  is  undefined. 

If  the  incident  neutron  energy  is  greater  than  Eg^,  two  cases  must  be  considered- 

p 

1.  If  Ej  is  sufficiently  large  that  the  probability  of  discrete  level  scattering 
is  zero,  then  only  continuum  level  scattering  occurs.  Let  ^BD  be  the 
energy  above  which  all  a^,  are  zero,  then  the  kernel  g(Ej  ,E')  can  be  rep¬ 
resented  by  a  continuous  probability  distribution  function,  g^(Ep,E') 
which  satisfies  the  following  condition. 

Ep 

gc(Ep,E')  dE'  =  1  (8) 

2.  If  Egp  <  Ep  ^  Egg,  then  both  continuum  and  discrete  level  scattering 
occur.  Both  the  continuum  and  discrete  kernels,  gg  and  a^’s  are  comput¬ 
ed,  satisfying  the  normalization  given  by  Eqs.  7  and  8.  Hence,  the  in¬ 
elastic  scattering  kernel  is 

gCEP.E  )  =  ‘'nn-,c(Ep)gc(Ef,E-)  ^  D(Ef)  gp(Ep,E-) 

'  ’  <'nn',c(Ep)  +»nn',D(Ep) 

where  o'nn'^c  ^nn',D  cross  sections  for  the  excitation  of 

the  continuum  and  of  the  discrete  levels,  respectively. 
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The  continuum  energy  spectrum  of  neutrons  emitted  per  inelastic  scattering  at 
energy  E  is  obtained  using  a  routine  called  ASPIC  which  computes  the  probability 
g(E,E')  according  to  a  statistical  theory.* 


g(^)(E,E') 


E^4^^  (E0p(^^(E-E^) 


E'aj3^(E')  p(^^(E-E')  dE' 


(10) 


where  cr^^^E)  is  the  cross  section  for  compound  nucleus  formation  and  P^^\e)  is 
the  level  density  in  the  nucleus  i  (i  =  1  refers  to  the  target  nucleus). 

The  following  form  is  assumed  for  the  level  density: 

p(^^(E)  =  0  for  E  <  Ep 

P^^^E)  =  po  for  eS^^<  E  <  eJ^^  (11) 

p^^^(E)=  Po  exp  [2/ a^^^E  “  2/a(^^Eo^’-^  ]  for  E  >  Ep. 


For  the  treatment  of  (n,2n)  reactions,  it  is  assumed  that  the  first  neutron  has  the 
distribution  (Eq.  10)  where  i  =  1.  After  the  emission  of  the  first  neutron,  the  tar¬ 
get  nucleus  is  left  at  excited  states  E"*’  with  a  population  density: 


R(E,E+)  =  g(i)(E,E-E+).  (12) 

R(E+)  has  a  maximum  at  E'*’-  E  -  /E/aS^  ['/E/a^  is  the  nuclear  temperature]. 

We  make  the  simplifying  assumption  that  all  secondary  neutrons  are  emitted  from 
the  states  E'*’  =  E  -  / E/a.^^ .  Applying  statistical  theory,  the  spectrum  of  second¬ 
ary  neutrons  per  (n,2n)  collision  is  given  by 

fE-  /z/Wi  -  E^,E']  (13) 

*See  terminology  and  definitions  at  the  end  of  this  section. 


14 


where  is  the  threshold  for  the  (n,2n)  reaction,  is  given  by  Eq.  10,  and  the 
superscript  (2)  refers  to  the  nucleus  left  after  the  (n,2n)  reaction  on  the  target. 

In  order  to  account  for  the  neutrons  leading  to  the  ground  state  of  this  final  nucle¬ 
us,  one  should  set  e[^^  ~  0  [choose  e|^^  =  0.01  J.  Hence,  the  composite  spectrum 
of  neutrons  emitted  per  (nn';  n,2n)  reaction  is  given  by 

g(E,E')  =  g(‘)(E,E')  +  g^^)(E,E')  (14) 

'^nn'  ,c^E^ 

where  included  into  GENDA. 

In  the  event  that  the  core  storage  requires  more  than  18,000  words  for  g(E,E'), 
the  code  determines  a  lowest  value  for  E'  (which  will  be  higher  than  the  lowest 
energy  of  the  problem)  so  that  g(E,E')  fits  in  storage.  The  portion  of  g(E,E') 
neglected  is  generally  negligible. 

The  above  description  of  the  calculation  of  the  probability  g(E,E')  is  not  adequate 
for  light  nuclei.  Special  treatments  for  lithium  and  beryllium  are  available  at 
the  user’s  option. 

The  GENDA  Li  option^  (see  Fig.  2)  .assumes  a  continuous  spectrum  for  Ej  ^ 

3.95  Mev.  For  E  below  3.95  Mev,  the  inelastic  scattering  arises  only  from  the 
excitation  of  the  0.477-Mev  level. 

The  continuous  probability  distribution  function,  gj,(E,E'),  is  taken  as  constant  with 
a  value  of 

^max  =  (1^^ 

and  stretches  from  E'  =  2  x  10'®  Mev  to  an  E,'  such  that  E'.  ^  E^„^  ?  E,' 

j  j+i  max  j 

where  Ej  and  Ej^j  are  adjacent  energies  in  the  integration  mesh.  It  should 


15 


6 


-  Lithium  g(E,  E')  computation 


This  Document  Contains 
Missing  Page/s  That  Are 
Unavailable  In  The 
Original  Document 


1^lc5unio~py 

buLMCr^ 

'H^fVC'ieh 


BEST 

AVAILABLE  COPY 


be  noted  that 


I  g^(E.E')  dE'  -  -t  (16) 

^max 

is  not  equal  to  1,  though  close  to  it.  Because  of  the  Li  nucleus  recoil,  gQ(E,E') 
is  obtained  as  follows: 


gj)(E,E') 


2  6  (E'~E| ) 


(17) 


where  for  each  E  the  value  of  Ej  is  selected  which  lies  closest  to 

Eq  =  0.781  E  -  0.420  Mev  (18) 

The  GENDA  Be  option^  (see  Fig.  3)  uses  revised  data  as  obtained  by  Goldstein 
and  Mechanic.^ 


Briefly, 

1.  For  eP  <  2  Mev: 

J 

g(EG,  E')  =  0 

2.  For  2  Mev  ^  eP  2.83  Mev: 

Of  the  two  types  of  neutrons  emitted,  ^(Ej  )  go  via  excitation  of  the 
2.43-Mev  level  of  Be®.  The  remaining  2  -  e  neutrons  are  emitted  with 
a  constant  distribution  stretching  from  zero  to  the  highest  energy  allowed 
(namely  E  -  1.7). 


gc(Ef,E')  = 


2  -  e(Ep) 


(19) 


3.  For  eP  >  2.83  Mev: 

The  contribution  from  the  discrete  level  must  be  added.  Thus,  the  over¬ 
all  distribution  g(Ep,E')  becomes 
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g(Ep,E')  =  g^,(Ep,EO  +  €{E^) 
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1 


10V0.8  E^-2.187 

- w — 

] 


2.5 


(20) 


€(Ej^  is  taken  from  the  data  compiled  by  Lustig  and  Kalos. 


Terminology  for  Section  2. 1 


^t 

a(‘),  eS^^ 

eS‘^ 

eS^^ 

a(0 


the  threshold  of  the  (n,2n)  reaction 
level  density  parameters  for  the  tar  get  nucleus 
level  density  parameters  for  the  nucleus 
energy  at  which  the  continuum  starts 
should  be  set  equal  to  0.01 

cross  section  for  ccmpou.nd  nucleus  formation  of  the  target 
nucleus 

cross  section  for  compound  nucleus  left  over  after  the  (n,2n) 
reaction. 


2.2  THE  PREPARATION  OF  INPUT  DATA  FOR  GENDA 

T^nless  otherwise  stated,  all  data  for  this  code  are  in  standard  FORTRAN  E  for¬ 
mat.  The  data  are  punched  in  six  fields  of  10  columns  each  (i.e.,  6E10.4).  For 
the  general  data  decks,  it  is  in  the  order  of  energy  followed  by  its  cross  section. 


Example 


Column 

1-10 

El 

11-20 

a(Ei) 

21-30 

Ez 

31-40 

a(E2) 

41-50 

Eg 

51-60 

a(E3) 

10 


20 


30 


40 


50 


60 


/ 

Et 

a(Ei) 

E2 

^(Ej) 

Ej 

a(Es) 

All  energies  must  span  the  range  from  to  (the  minimum  energy) 

with  the  exception  of  the  discrete  inelastic  scattering  cross  section  (cTnn'^L). 

The  data  control  information  is  in  integer  format,  with  the  exception  of  the 
true  mass  number  which  is  in  F  format. 

2.3  GENDA  PROBLEM  INPUT 

2. 3. 1  Title  Card 

Any  80  Hollerith  characters  (e.g.,  user’s  name,  element  name,  date,  etc.) 

2.3.2  Problem  Control  Card 


Column  Item 

1-3  Number  of  elements  (NEL) 

4-7  Number  of  final  energies  (JMAX)  (<450) 

8  lEN  =  0  -  compute  final  energy  table  with  fixed  lethargy 
lEN  =  1  -  final  energy  table  is  input 

9  IPHI  =  0  -  weighting  function  0^(E)  is  not  present 
IPHI  =  1  -  weighting  function  4>(.{E)  is  present 

10  IDELP  =  0  -  table  of  B^E,-)  is  not  present 
IDELP  =  1  -  table  of  B^E^)  is  present 

Note: 

IDELP  must  =  0  if  IPHI  =  0. 
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2.3.3  Energy  Input 


For  lEN  =  0 

Column 

1-10 

11-20 


Item 


Emax  “  ^^Aimum  energy  (Mev) 
H,  fixed  lethargy  step 


1 

JMAX 


1„  Emm 

ejmax 


Floating 

Floating 


For  lEN  =  1 

Energy  table  eP,  j  =  1,  2,  ...,  JMAX;  Ep  >  eP 
Format  (6E10.4) 


2.3.4  Weighting  Function  0p(E)  (IPHl  -  1) 


Identification  Card 


Column 


Item 


1-6  Any  six  Hollerith  characters  for  identification 

8  =  1  -  linear 

-  2  -  linear  on  log-log  paper 
9-12  Number  of  (p^{E)  tabulated  (^300) 

E*^,  Data  Deck 


Note: 

The  must  be  included  if  CAR  type  group  averaging  is  to  be  used. 


2.3.5  Weighting  Function  B^(E^)  (IPHI  =  1  and  IDELP  -  1) 


Identification  Card 


Column  Item 

1-6  Any  six  Hollerith  characters 

7-8  Number  of  B-j.  tabulated  (<99  entries) 


E^  Bt(Et)  Data  Deck 
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2.4  GENDA  ELEMENT  INPUT 


2.4.1  Element  Control  Card 


Column 


Item 


1-6 

7-13 

14-23 

DDDDZZZEEE 
14  23 


24-50 


24 

25 

26 

27,  28 
29 

30,  31 
32 

33,34 

35 

36 

37 
3d 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 


Elemojr.  name,  Hollerith  format 
True  mass  (F  format) 

Integral  mass  (numerical  identification)  integer. 

D  -  field  for  data 
Z  -  field  for  atomic  number 
E  -  field  for  element  identification 
E  =  000  if  natural  element 
Control  for  the  following 
Data  =  1  -  data  are  present 
Data  =  0  -  data  are  not  present 
Blank 

Total  cross  section,  Orp 

Elastic  scattering  cross  section, 

The  order  of  expansion  (T^  <  20,  f, 

elastic  scattering  is  anisotropic 

Excitation  of  discrete  levels,  t 

T  nn  jJj 

The  order  of  expansion  (T^)  if  the  excitation  of  discrete  levels 
is  anisotropic 

Excitation  of  the  continuum,  ^  (includes  (Tn,2n) 

The  order  of  expansion  (T^)  if  the  excitation  of  the  continuum 
is  anisotropic 


P’ 


p  =  1,  2,  ...,  tN)  if  the 


Cross  section  for  the  emission  of  two  neutrons  o. 


n,2n 


Fission  cross  section,  Of 
Absorption  cross  section, 

Transport  cross  section,  cr^j. 

Cross  section  for  the  emission  of  an  o;  particle 
Cross  section  for  the  emission  of  a  proton,  a 
Cross  section  for  the  emission  of  Hj,  ^ 

Cross  section  for  the  emission  of  Hj, 

Radiative  capture  cross  section,  a, 

Fission  y  spectra 
Absorption  y  spectra 
Nonelastic  y  spectra 
t'(E)  data 
77(E)  data 

Edit  synthesized  fp  data 
Thermal  data 


n,y 


2.4.2  Total  Cross  Section, 


Identification  Card 


Column 

Item 

1-6 

Any  six  Hollerith  characters 

7 

=  1  —  linear 

=  2  -  linear  on  log-log  paper 

10 

=  0  -  <7(E^)  found  by  interpolation 
=  1  —  (^Ej  )  found  by  CAR  type  averaging 
Number  of  a(E)  tabulated  (<1100) 

11-14 

15 

=  0  -  no  resonance  parameters 
=  1  -  resonance  parameters,  for 

E,  arp(E)  Data  Deck 

Thermal  Value  of  (Format  El 0.4) 
Identification  Card  for 


Column 

Item 

1-6 

Any  six  Hollerith  characters 

7-8 

Number  of  Aj^  tabulated  (^99 

Ap(Ei^)  Data  Deck 


2.4.3  Elastic  Scattering  Cross  Section,  (?„ 


Identific^’*'.on  Card 


Column 

Item 

1-6 

Any  six  Hollerith  characters 

7 

=  1  -  linear 

=  2  -  linear  on  log-log  paper 

10 

=  0  -  cr(Ej^  found  by  interpolation 
=  1  -  CT  (Ep)  found  by  CAR  type  averaging 
Number  ol  a(E)  tabulated  (<1100) 

11-14 

15 

=  0  -  no  resonance  parameters 
=  1  -  resonance  parameters,  A^^  for 

E,  o^jE)  Data  Deck 


Thermal  Value  for  (Format  E10.4) 
Identification  Card  for 


Column  Item 

1-6  Any  six  Hollerith  characters 

7-8  Number  of  A^,  tabulated  (<99) 

E;,,  A„(Et,)  Data  Deck 

2.4.4  Angular  Distribution  Data,  f^ 

Anisotropic  elastic  scattering  for  each  fp,  p  =  1,  2,  ...,  (T^  ^  20). 

Identification  Card  (one  for  each  fp  data  deck) 

Item 

Any  six  Hollerith  characters 
=  1  -  linear 

=  2  —  linear  on  log -log  paper 
=  0  -  fp(Ep)  found  by  interpolation 
~  1  ~  Ip(Ep)  found  by  CAR  type  averaging 
Number  of  fp(E)  tabulated  (-*?  1100) 

2.4.5  Inelastic  Scattering  Data,  Discrete,  On,/  I. 

Inelastic  Scattering  Data,  Discrete,  Control  Card 

Item 

Ebd»  energy  at  which  the  discrete  spectrum  begins.  This 
value  can  be  less  than  E^^  (e.g.,  Ej^^  =  18.01739  Mev, 
EgD  "  7.93  Mev) 

NU,  number  of  level  parameter-,  L  <"20) 


Column 

1-10 

11-13 


Column 

1-6 

7 

10 

11-14 
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Excitations  of  Discrete  Levels,  L 
Identification  Card 


Column 

Item 

1-6 

Any  six  Hollerith  characters 

7 

=  1  -  linear 

=  2  -  linear  on  log-log  paper 

10 

=  0  -  a  (i^)  found  by  interpolation 
=  1  -  a(Ej  )  found  by  CAR  type  averaging 

11-14 

Number  oi  a  (E)  tabulated  (^1100) 

15 

=  0  -  no  resonance  parameters 
=  1  -  resonance  parameters,  Aj^,  for 

^BD  max 

Thermal  Value  of 

%n',L  (Format  E10.4) 

Identification  Card  for  Aj^ 

Column 

Item 

1-6 

Any  six  Hollerith  characters 

7-8 

Number  of  A^,  tabulated  (<99) 

E^,  Aj/(Ej^)  Data  Deck 
^  Deck 

u  =  I,  2y  NU  Ej  <  E2  <  <  Ej^  (format  6E10.4) 

Level  Parameters  ^ 

Identification  Card  (one  for  each  deck  of  1^=  1,2, 

Column  Item 

1-6  Any  six  Hollerith  characters 

7  =  1  -  linear 

=  2  -  linear  on  log -log  paper 


Column 

Item 

10 

=  0  -  L(Ep)  found  by  interpolation 

=  1  -  (Ep)  found  by  CAR  type  averaging 

Number  of  ^u^nn',L  tabulated  (<1100) 

11-14 

Deck  from  EgD  to 

Thermal  Value  of  (Format  E10.4) 

Identification  Card  for  Aj^ 

Column  Item 

1-6  Any  six  Hollerith  characters 

7-8  Number  of  Aj,  tabulated  (<99) 

Ep,,  Aj;(Ej^)  Data  Deck 

Anisotropic^Oj^j^^^T 

The  routines  to  handle  this  cross  section  have  not  been  programmed. 
2.4.6  Excitation  of  the  Continuum,  p  (Includes  gn^m) 


Inelastic  Scattering  Ekita,  Continuous,  Control  Card 


Column  Item 

1-10  Egc 

13  IGOP:  Option  to  commute  G(E,E') 

IGOP  =  1  -  ASPIC 
IGOP  =  2  -  lithium  special  code 
IGOP  =  3  -  beryllium  special  code 
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Excitation  of  the  Continuum, 


a 


nn^c 


Identification  Card 

Item 

Any  six  Hollerith  characters 
=  1  -  linear 

=  2  -  linear  on  log-log  paper 
=  0  -  a(Ep)  found  by  interpolation 
=  1  -  or(Ep)  found  by  CAR  type  group  averaging 
Number  of  cr(E)  tabulated  (^1100) 

=  0  -  no  resonance  parameters 
=  1  -  resonance  parameters,  for  <^nn',c 

^nn',c^^^  (^MAX  ""  ^MIN^ 

Thermal  Value  of  a  /  (Format  E10.4) 

till  J  ^ 

Identification  Card  for  Aj^ 

Column  Item 

1-6  Any  six  Hollerith  characters 

7-8  Number  of  tabulated  (<99) 

E^,,  ApiE^,)  Data  Deck 


Column 

1-6 

7 

10 

11-14 

15 


IGOP  =  1 


Parameters  for  the  Temperature  Model  of  the  Nucleus 


Column 


Item 


1-10 

11-20 

21-30 


31-40 

41-50 

51-60 

61-70 


Ei,  the  threshold  energy  of  the  (r.,2n)  reaction 

a^^^  -  level  density  parameter  for  the  target  nucleus 

a^^^  -  level  density  parameter  for  the  nucleus  after  an 

!n,2n)  reaction 

level  density  parameter  for  the  target  nucleus 
level  density  parameter  for  the  nucleus  after  an 
fn,2n)  reaction 
e(^^  -  first  excited  state 
Ej^^  -  should  be  set  equal  to  0.0 
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Cross  Section  for  the  Emission  of  Two  Neutrons  a 

’  n,2n 

1.  Identification  Card 


Column 


Item 


Any  six  Hollerith  characters 
=  1  -  linear 

on  log-log  paper 

found  by  interpolation 

found  by  CAR  tjrpe  group  averaging 

o<E)  tabulated  (^1100) 

5  =  0  -  no  resonance  parameters 

=  1  -•  resonance  parameters,  A„  for  a 

^  n,2n 

2.  E,  ,„(E)  Data  Deck 


=  2  -  linear 
10  =  0  -  a(E^) 

=  1  -  a(Ep) 
11-14  Number  of  ^ 


Value  of  (Format  E10.4) 
4.  Identification  Card  for 


Column 

1-6  Any  six  Hollerith  characters 

7'8  Number  of  A^  tabulated  (^9) 


5.  Et^,At/(Et/)  Data  Deck 


Cross  Section  for 

Compound  Nucleus  Formation  of  the  Target  Nuc: 

1.  Identification  Card 

Column 

Item 

1-6 

Any  six  Hollerith  characters 

7 

=  1  -  ’  inear 

10 

=  2  -  linear  on  log -log  paper 
~  0  —  a(Ei  )  found  by  interpolation 

1  —  <7(Ej  )  found  by  CAR  type  group  averaging 
Number  of  o(E)  tabulated  (^1100) 

11-14 

2.  E,  a<^)(E)  Data  Deck 


Cross  Section  for  Compound  Nucleus  Left  Over  After  the  (n,2n)  Reaction, 

1.  Identification  Card 

Column  Item 

I- 6  Any  six  Hollerith  characters 

7  =  1  -  linear 

=  2  -  linear  on  log -log  paper 
10  =  0  -  CT(Ej  )  found  by  interpolation 

=  1-0  (eP)  found  by  CAR  type  group  averaging 

II- 14  Number  ol  a(E)  tabulated  (^1100) 

2.  E,  Data  Deck 

IGOP  =  2  (  Lithium  Special  Code) 


No  input  required. 


IGOP  =  3  (Beryllium  Special  Code) 


Identification  Card 

Item 

Any  six  Hollerith  characters 
=  1  -  linear 

=  2  -  linear  on  log -log  paper 
=  0  -  €  (^)  found  by  interpolation 
=  1  -  e(Ep)  found  by  CAR  type  group  averaging 
Number  of  e(E)  tabulated  (<1100) 

E,  e(E)  Data  Deck 


Column 

1-6 

7 

10 

11-14 


Anisotropic  o  , 

nil  ^ 

The  routines  to  handle  this  cross  section  have  not  been  programmed. 
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2.4.7  Cross  Section  for  the  Emission  of  Two  Neutrons,  a  . 
_ *  n,2n 

(If  <7n^2n  supplied  for  the  ASPIC  routine  ^  not  input  again  and  ^not  flag 

^n,2n  element  control  card.) 


Identification  Card 


Column  Item 


1-6 

Any  six  Hollerith  characters 

7 

=  1  -  linear 

=  2  -  linear  on  log-log  paper 

10 

=  0  -  found  by  interpolation 

=  1  -  a  (eP)  found  by  CAR  type  group  averaging 

11-14 

Number  oi  o(E)  tabulated  (^1100) 

15 

=  0  -  no  resonance  parameters 
=  1  -  resonance  parameters,  for 

E>gn.2n(E)  Data  Deck 

Thermal  Value  of  an.2n  (Format  E10.4) 

Identification  Card  for 

Column  Item 

1“6  Any  six  Hollerith  characters 

7-8  Number  of  Aj/  tabulated  (^99) 

Ei;,  At/(Et/)  Data  Deck 

For  the  following  eight  cross  sections  (see  Section  2.4.8)  the  following  general 
format  is  to  be  used  in  the  order  of  the  cross  sections  listed. 

Identification  Card 

Column  Item 

1-6  Any  six  Hollerith  characters 

7  =  1  -  linear 

=  2  -  linear  on  log -log  paper 
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10  =  0  -  cr(Ej^  found  by  interpolation 

=  1-0  (Ej  )  found  by  CAR  type  group  averaging 
11-14  Number  of  a(E)  tabulated  (<1100) 

15  =  0  -  no  resonance  parameters 

=  1  -  resonance  parameters,  Aj^  for  a 


E,  a<E)  Data  Deck 

Thermal  Value  of  a  (Format  E10.4) 

Identification  Card  for 

Column  Item 

1-6  Any  six  Hollerith  characters 

7-8  Number  of  A jy  tabulated  (^99) 

Ei^,  At^(Ej^)  Data  Deck 

2.4.8  Other  Cross  Sections 

1.  Fission  Cross  Section,  of 

2.  Absorption  Cross  Section,  o^^j^g 

3.  Transport  Cross  Section,  o^^ 

4.  Cross  Section  for  the  Emission  of  an  a  Particle, 

5.  Cross  Section  for  the  Emission  of  a  Proton,  Op 

6.  Cross  Section  for  the  Emission  of  H2, 

7.  Cross  Section  for  the  Emission  of  H3, 

8.  Radiative  Capture  Cross  Section, 

The  routines  to  handle  fission  y  spectra,  absorption  y  spectra,  nonelastic  y  spec¬ 
tra,  i'(E)  data,  and  rj(E)  data  have  not  been  programmed. 

If  the  program  control  card  specifies  more  than  one  element,  repeat  Section  2.4 
of  the  above  input  preparation  for  each  element. 
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2.5  GENOA  OUTPUT 


The  GENOA  output  consists  of  two  edits.  The  first  edit  tabulates  the  input  data 
in  the  same  order  as  they  appear  in  the  input  cards.  The  second  edit  tabulates 
the  cross-section  data  at  the  output  mesh  energies,  as  described  in  Section  2.1, 
and  arranges  them  in  an  easily  readable  manner. 
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3.  THE  GENERATION  DATA  ROUTINE  PROCESSOR  (GENPRO) 


The  cross-section  data  generated  by  GENDA  may  be  used  in  various  ways.  How¬ 
ever,  in  Monte  Carlo  type  calculations,  what  has  to  be  determined  once  a  collision 
occurs, is 

1.  The  type  of  event  which  takes  place 

2.  The  energy  and  direction  angle  of  the  emerging  particle,  if  any. 

The  GENPRO  program  processes  the  GENDA  output  to  tabulate  probabilities  pecul¬ 
iar  to  a  nuclide  at  a  prescribed  energy  mesh  as  follows. 

1.  It  computes  the  probability  of  elastic  scattering,  inelastic  scattering,  and 
absorption. 

2.  In  case  of  anisotropic  elastic  scattering,  it  tabulates  the  cosines  of  angles 
between  which  scattering  is  equi-probable. 

3.  In  case  of  inelastic  scattering,  two  possibilities  exist. 

a.  If  the  spectrum  has  a  continuous  part,  it  computes  energies  to  which 
the  particle  scatters  with  equal  probability. 

b.  If  the  spectrum  is  representable  by  v  distinct  levels,  it  tabulates  the 
probability  that  the  particle  will  scatter  to  the  level  E^. 

It  is  important  to  note  at  this  point  that  the  maximum  number  of  output  energies  is 
different  for  GENDA,  GENPRO,  and  the  subsequent  Monte  Carlo  programs.  In¬ 
deed,  because  of  the  need  of  fine  tabulation  of  particular  cross  sections,  GENDA 
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was  expanded  to  handle  450  output  energy  points.  On  the  other  hanc,  GENPRO  can 
only  handle  401  energy  points  due  to  core  storage  limitations.  The  U>tc  Monte 
Carlo  programs  are  further  restricted  to  a  maximum  energy  mesh  of  81  points. 
Therefore,  in  any  Monte  Carlo  type  calculations,  input  to  GENDA  must  be  consist¬ 
ent  with  the  above  restrictions. 

Also,  the  output  energy  mesh  for  GENDA  is  variable.  However,  if  it  is  used  in 
connection  with  GENPRO,  a  fixed  lethargy  interval  must  be  used. 

3.1  TECHNICAL  DESCRIPTION 

The  cross-section  data  for  each  element  provided  by  GENDA  are  written  on  a 
magnetic  tape  called  the  Genda  Output  Tape  (GOT). 

GENPRO  uses  GOT  to  create  an  Element  Data  Tape  (EDT)  which  contains  the 
following  information: 

1.  The  number  of  energy  points,  J 

2.  The  set  of  final  energies  in  increasing  order  E?  (using  the  given  constant 
lethargy  step). 

Then,  for  each  element, 

1.  The  total  cross  section  at  each  final  energy,  ffrj.(Ep),  is  transmitted  to  the 
EDT 

2.  The  probabilities  of  elastic  and  inelastic  scattering  and  absorption  are 
computed. 


(21) 


(23) 


Pa(E°)  =  1  -  Ps(eG)  -  P^,(Ef) 

Pg  and  tables  are  then  transmitted  to  the  EDT. 

It  should  be  noted  that  for  neutron  absorption,  P;^  contains  both  fission 
and  capture  probabilities. 

3.  In  the  event  of  elastic  scattering,  the  following  calculation  is  performed 
to  decide  whether  or  not  the  scattering  is  isotropic. 

The  angular  distribution  data  supplied  by  GOT  in  the  form  of  Legendre 
expansion  coefficients  ff 's  is  tested  for  each  energy  point.  Let 

L 

a(w=l)  =  Y,  (2^  +1)  fjp(Ep) 

f=0 

L 

a(w  =  -l)=  2  (2f  +  l)fjg(Ep)  (-1)^ 

11=0 

a(a)=0)=  Y,  (2f +  1)  ff(Ep)  Pjp(O) 
even  f 's 

where  the  notation  used  is  given  in  the  GENDA  section. 

The  largest  and  smallest  of  these  are  selected  and  referred  to  as  M  and 
m,  respectively. 

Then  if 


(24) 


(25) 


(26) 


M  -  m 
0.5(M+m) 


^  0.05, 


(27) 


the  angular  distribution  is  taken  as  isotropic  and  no  further  tabulated 
data  are  required.  (See  Section  7.) 
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otherwise,  the  scattering  is  treated  as  anisotropic  and  GENPRO  proceeds 
to  compute  the  aiigles  6  through  which  it  is  equally  probable  to  scatter. 


Since  the  differential  scattering  cross  section  is  given  by 


<'(e”  m)  = 


gn(Ep) 

An 


L 

2  (2«  +  l)f,(Ep)PjM 
«=0 


where  p  =  cos  9  varies  from  -1  to  +1,  GENPRO  computes 


(28) 


J  •/-I 


(29) 


which  gives  the  p  probability  distribution.  The  integration  of  Eq.  28 
yields 


n  ^  1 

F(u,Ef)=  2  5  (Ep)  {Pf+,W  -  P{-,(m)1  (30) 

«=0 

with  the  convention  that  P_j  =  -1. 

F(ji,Ej^  is  computed  for  -l^/i<+l  in  steps  of  0.01.  Then  GENPRO  locates 
such  that 

F(Mn)  for  n  =  0,  1,  2  ...,  N  (31) 

N  is  usually  10  for  the  Monte  Carlo  programs. 

For  statistical  considerations,  it  is  desirable  to  use  a  new  variable  y 
which  is  related  to  p  by 


H  =  cos  0  =  1  -  2x 

where  now  x  varies  from  0  to  1. 

Hence, 

=  >.2-  N 

is  tabulated  and  transmitted  to  the  EDT  as  Table  CHI. 


(32) 


(33) 


4.  In  the  event  of  continuum  inelastic  scattering  (which  is  determined  by  the 
threshold  energy  supplied  by  GOT),  the  probability  that  a  neutron  enter- 

Q 

ing  collision  at  energy  will  emerge  from  an  inelastic  scattering  colli¬ 
sion  at  energy  E'  has  been  computed  by  GENOA  and  entered  into  GOT  in 
the  form  of  gj.(E9,E')  with  Ep  >  Egp. 


To  obtain  the  energies  to  which  it  is  equally  probable  to  scatter,  GENPRO 
calculates 


XE 

gc(Ep,E0dE',  E<Ep 


where  the  integral  is  evaluated  by  the  trapezoidal  rule. 

Note  that  if  one  considers  the  spectrum  of  neutrons  emitted  by  (nn')  re¬ 
actions,  the  integral  (Eq.  34)  is  normalized  to  1. 

Indeed, 

G(Ep,E9)  =  1  (35) 

^  J 

However,  GENOA  considers  the  composite  spectrum  of  neutrons  emitted 

G  G 

by  (nn';  n,2n)  reactions  and  hence,  G(Ej  ,Er)  may  be  different  from  1 
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and  GENPRO  proceeds  to  compute: 


«(eGe)  = 


G(Ef,E) 

G(EP,Ep) 


Then  it  locates  Ej^  such  that 
«(EGEn)=  ^forn  =  0,  1,  2  N 


(36) 


(37) 


N  is  usually  10  for  the  Monte  Carlo  programs.  The  Ej^  are  tabulated  and 
transmitted  to  the  EDT  as  Table  ENN. 


1.  In  the  event  of  discrete  level  inelastic  scattering,  GOT  supplies  the  aj,*s. 
The  a|,*s  are  tabulated  and  transmitted  to  the  EDT  as  Table  PLEV  (t*,]). 
Tf>e  corresponding  values  of  the  excitation  level  energies  are  tabulated 
am  transmitted  to  the  EDT  as  Table  ELEV  (i^.j). 


6.  If  both  continuum  and  discrete  level  inelastic  scattering  occur,  i.e., 

Egc  <  Ep  <  ^BD*  supplies  information  on  the  kernels  gc(Ep,E')  and 
g£)(Ep,E0  and  the  continuum  and  discrete  level  inelastic  scattering  cross 
sections,  ^nn',D^®'P^* 


To  obtain  the  energies  to  which  it  is  equally  probable  to  scatter,  the 
GENPRO  program  calculates  the  probability  G(Ep,E)  that  a  neutron  at 
Ep  scatters  at  energy  E  as  follows: 


rE?  rEf 

W,c(Ep)  j  gc(Ep,E')dE'-.a^,  D(Ep)  J  gD(Ep,E0dE 


G(Ep,E)  = 


W,c(Ef)  -  o^^^D(Ef) 


(38) 
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where  G(eP,0)  =  0  and  G(Ep,Ep)  =  1,  and  the  probability  distribution 
functions  and  gj)  have  been  discussed  in  the  GENDA  section. 

Then,  GEIJPRO  proceeds  to  locate  the  energies,  Ej^,  to  which  it  is  equally 
probable  to  scatter  as  if  the  level  spectrum  was  continuous  only.  The 
method  used  is  similar  to  that  of  determining  the  cosines  of  the  angles 
through  wkich  it  is  equally  probable  to  scatter.  The  En  are  tabulated  and 
transmitted  to  the  EDT  as  Table  ENN. 

3.2  GENPRO  INPOT 

The  input  to  GENPRO  consists  of  the  GENDA  Output  Tape  (GOT)  (Logical  tape  9) 
and  one  card  denoting;  the  number  of  elements  on  the  GOT  in  the  format  110. 

3.3  FORMAT  OF  THE  EDT 

A  description  of  the;  format  of  the  EDT  follows.  Since  the  EDT  is  made  up  of 
card  images,  the  dejscription  is  in  terms  of  card  images. 

Item  No.  oi' 

No.  Entries  Card  Description  Format 

1  1  NENERG  110 

NENERG  is  the  number  of  energies  in  the 
energy  table. 

2  NENERG  ETABLE(l),  ETABLE(2),  ETABLE(3),  5E14.5 

ETABLE(4),  ETABLE(5),  ... 

ETABLE  (NENERG). 

ETABLE  is  the  energy  mesh  for  all  ele¬ 
ments.  ETABLE(l)  is  the  lowest  energy 
and  ETABLE  (NENERG)  is  the  highest 
energy. 

3  3  AWT.  lAWT,  J  El 6. 8,  216 

AWT  is  the  floating  point  atomic  weight. 
lAWT  is  the  fixed  point  atomic  weight. 

J  is  an  end  of  data  flag.  J^^O  means  this 
card  is  the  last  card  of  the  EDT. 
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Item 

No, 

No.  of 
Entries 

Card  Description 

Format 

4 

5 

OF,  NOW,  NFI,  ISI,  IPIN 

OF^tO  means  this  is  last  card  of  element. 
NOW  is  the  number  of  words  in  the  next 
group  of  data.  NFI  is  a  file  number  used 
for  identification  only.  ISI^tO  means  the 
next  group  of  data  is  the  scatter  index. 
IPIN^O  means  this  element  has  no  inelas¬ 
tic  scattering. 

E16.8,  416 

5 

NOW 

SIGMA (1),  SIGMA (2),  SIGMA(3),  SIGMA(4), 
etc. 

SIGMA  is  the  microscopic  total  cross  sec¬ 
tion.  Note  that  NOW  must  equal  NENERG 
and  that  OF,  ISI,  IPIN  must  all  be  zero. 

5E14.5 

6 

5 

Repeat  item  4. 

S16.8,  416 

7 

NOW 

PSCAT(l),  PSCAT(2),  PSCAT(3),  PSCAT(4), 
PSCAT(5),  ...,  PSCAT (NENERG). 

PSCAT  is  the  probability  of  elastic  scatter¬ 
ing.  PSCAT(l)  corresponds  to  ETABLE(l) 
and  PSCAT  (NENERG)  corresponds  to 
ETABLE  (NENERG). 

5E14.5 

8 

5 

Repeat  item  4. 

E16.8,  416 

9 

NOW 

PABS(l),  PABS(2),  PABS(3),  PABS(4), 
PABS(5),  ...,  PABS(NENERG). 

5E14.5 

10 

5 

OF,  NOW,  NFI,  ISI,  IPIN 

ISI  must  equal  1  because  the  next  data 
group  is  the  scatter  index.  If  this  ele¬ 
ment  has  no  inelastic  scatterir^,  then 

IPIN  roust  equal  1. 

El 6. 8,  416 

11 

NENERG 

IP,  ID,  IT,  lA 

There  are  NENERG  cards  in  this  data 
group.  Each  card  has  the  quantities  IP, 

ID,  IT,  lA  defined  as  follows: 

4110 

IP  =  1  for  inelastic  discrete  scattering 
=  2  for  inelastic  continuum  scattering 

ID  =  the  location  of  an  ENN  or  PLEV 

table  relative  to  the  first  word  of  the 
SIGMA  table 
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Item  No.  of 

No.  Entries  Card  Description  Format 

11  NENERG  IT  =  1  for  isotropic  scattering  in  center  of  4110 

(Cont.)  mass 

=  2  for  scattering  in  hydrogen 
=  3  for  anisotropic  scattering  in  center 
of  mass 

=  4  isotropic  scattering  in  the  lab  system 
ID  =  location  of  a  x  Table  if  IT  =  3. 

12  5  Repeat  item  4  with  NOW  =  11.  E16.8,  416 

13  11  CHI(l),  CHI(2),  CHI(3),  CHI(4),  CHI(5),  etc.  5E14.5 

CHI  is  a  table  of  11  entries. 

CHI(l)  =  1.0,  CHI(ll)  =  0.  Items  12,  13 
are  repeated  for  each  energy  at  which 
“IT”  of  item  11  is  3. 

14  5  Repeat  item  4  with  NOW  =  11.  E16.8,  416 

15  11  ENN(l),  ENN(2),  ENN(3),  ENN(4),  ENN(5),  5E14.5 

etc. 

The  ENN  table  is  used  to  determine  the  en¬ 
ergy  after  scattering  for  an  inelastic  con¬ 
tinuum  interaction.  Items  14  and  15  are  re¬ 
peated  for  each  energy  at  which  ‘TP”  of 
item  11  is  2.  If  the  element  has  no  contin¬ 
uum  scattering  items,  14  and  15  are  omitted. 

16  5  Repeat  item  4  with  NOW  =  number  of  excita-  El 6. 8,  416 

tion  levels  for  discrete  scattering. 

17  NOW  ELEV(l),  ELEV(2),  ELEV(3),  ELEV(4),  5E14.5 

ELEV(5),  ...,  ELEV(NOW). 

The  ELEV  tabic  is  a  list  of  possible  excita- 
ticMi  levels;  each  entry  corresponds  to  a 
probability  in  a  PLEV  table. 

18  5  Repeat  item  16.  E16.8,  416 

19  NOW  PLEV(l),  PLEV(2),  PLEV(3),  PLEV(4),  5E14.5 

PLEV(5),  ...,  PLEV(NOW). 

The  PLEV  tables  are  tables  of  probabil¬ 
ities  of  scattering  from  the  current  en¬ 
ergy  to  the  excitation  levels.  Items  18 
and  19  are  repeated  for  bach  energy  at 
which  inelastic  discrete  scattering  oc¬ 
curs  (IP=1). 
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Item 

No. 

20 

21 


No.  of 

Entries  Card  Description 

5  OF,  NOW,  ISI,  IPIN 

OF  =  1.0.  NOW,  ISI,  IPIN  =  0.  Item  20 
is  the  last  card  of  an  element. 

3  AWT,  lATWT,  J 

If  there  are  no  more  elements  AWT, 
LATWT  =  0.  J  =  1.  If  there  are  more 
elements  repeat  from  item  3. 


Format 
El 6. 8,  416 
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4.  the  data  organization  program  (datorg) 


The  GENOA  and  GENPRO  programs  are  used  to  generate  an  element  data  tape 
(EDT)  which  contains  microscopic  cross  sections  and  probability  tables  for  all 
nuclides  of  interest.  This  EDT  can  be  considered  as  a  library  tape  for  any  sub- 
sequ^.  ^  Monte  Carlo  problems. 


For  a  specific  problem,  however,  the  system  is  divided  into  regions  which  contain 
given  nuclide  compositions.  Each  composition  may  be  described  by  a  set  of  nu¬ 
clides  and  their  respective  concentrations.  When  a  particle  enters  collision  at 
some  point  within  a  region  of  given  composition,  the  Monte  Carlo  program  has  to 
select  the  nuclide  in  the  composition  with  which  the  particle  collides. 

The  probability  that  a  particle  with  energy  E  will  interact  with  nuclide  k  of  the 
composition  R  is  given  by 


(39) 


where  j^(E)  is  the  total  macroscopic  cross  section  of  the  composition,  and 
Cr  is  the  concentration  of  nuclide  k  in  region  R. 


Ck,R  =  ^NoVk^R 


(40) 
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where  =  density  of  the  nuclide  k 
Aj^  -  its  atomic  weight 

No  10^^  =  0.6023  X  10^^,  the  Avogadro^s  number 
vjj  R  =  volume  fraction  of  the  nuclide. 

4.1  TECHNICAL  DESCRIPTION 

The  DATORG  program  searches  and  reads  from  the  EDT  all  information  concern¬ 
ing  the  nuclides  which  appear  in  the  various  coii  positions  specified  in  the  input. 

It  computes  the  total  macroscopic  cross  sections  of  each  composition  at  each  final 
energy,  using  the  concentration  data  supplied  as  input  to  the  program. 

Finally,  it  stores  in  the  COMMON  portion  of  core  all  data  concerning  the  nuclides 
and  the  composition  of  the  problem  in  a  prescribed  order. 

Hence,  the  role  of  DATORG  is  above  all  the  organization  of  the  cross  section  and 
probability  data  read  from  the  EDT.  In  addition,  it  computes  and  stores  the 
Pq'^R(Ep)^s. 


E 

k»sinR 


(41) 


The  output  energy  mesh  for  DATORG  is,  of  course,  given  by  the  fixed  lethargy 
interval  used  to  generate  the  EDT. 

4.2  DATORG  INPUT  FORM 

A  DATORG  input  form  is  shown  on  the  following  page. 
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This  Doctiment  Contains 
Missing  Page/s  That  Are 
Unavailable  In  The 
Original  Document 


T^lo-n  Ic,-  py 
iincd^  hai^ 


BEST 

AVAIUBLE  COPY 


numbered  seijuentUUy. 


4.3  DESCRIPTION  OF  THE  OUTPUT  DATA 

The  data  stored  by  DATORG  contains  the  following  items  in  their  order  of  listing. 


No.  of 

Symbol 

Locations 

Definition 

Counters 

NORG 

1 

Counter  denoting  number  of  words  in  the  array 
XSECT. 

NELEM 

1 

Counter  denoting  number  of  distinct  nuclides 
(«33). 

NR 

1 

Counter  denoting  number  of  compositions  (=s20). 

NENERG 

1 

Counter  denoting  number  of  energies  (<100). 

JK 

1 

Counter  denoting  the  sum  of  all  nuclides  in 
every  composition. 

Tables 

NEL 

NR 

Number  of  nuclides  for  each  composition  (<5 
for  any  composition). 

NEVAT 

NELEM 

Integral  value  of  atomic  weight  of  each  nuclide 

desired.  Elements  arranged  in  the  order 
read  from  EDT  (i.e.,  if  oxygen  =  16.0,  hydro¬ 
gen  =  1.0,  Fe  =  55.85  are  desired  and  they  ap¬ 
pear  on  the  EDT  as  16,  55.85,  1.0,  then 
NEL(l)  =  16,  NEL(2)  =  55,  NEL(3)  =  1.  Fur¬ 
thermore,  oxygen  will  be  referred  to  as 
element  1;  Fe,  as  element  2,  etc.). 


EVAT 

NELEM 

Corresponding  to  each  element  in  NEVAT  the 

NR 

actual  atomic  weight  is  entered  in  EVAT. 

lATWT 

£  NELi 

lATWT  contains  for  each  composition  in  the 

i=l 

order  of  composition  number  the  element 
number  (see  NEVAT)  of  each  element  of  the 
composition  in  the  order  specified  by  the  in¬ 

NR 

put. 

CONCT 

£  NELi 

Corresponding  to  each  entry  in  lATWT  the 

i=l 

concentration  of  the  element  is  placed  in 
CONCT. 

KPHYS 

NR 

There  is  one  entry  for  each  composition  and 
each  entry  contains  the  N^  w  here  Nj  =  Nj  _  j  + 

NELj.  J  +  1,  NELo  =  No  =  c. 
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No.  of 

Symbol 

Locations 

Definition 

XSECT 

(Variable 

dimension) 

Contains  the  cross  section  and  the  probabilities 
for  each  prescribed  element.  They  are  given 
in  the  order  in  which  they  are  obtained  in 
GENPRO  (and  appear  in  the  EDT).  Also  con¬ 
tains  macroscopic  cross  sections  for  each 
composition  in  the  order  of  composition  num¬ 
ber. 

NSECT 

NELEM 

Contains  the  origin  relative  to  start  of  XSECT 
for  each  element  (in  the  order  specified  by 
NEVAT). 

MUTORG 

NR 

Contains  origin  of  location  of  the  tables  of 
total  cross  sections,  in  XSECT  which  are 

written  in  the  order  of  composition  number. 
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5.  THE  GENERATION  OF  SOURCE  PARTICLES  FOR 
MONTE  CARLO  PROBLEMS 


The  Monte  Carlo  Game  played  in  the  solution  of  the  neutron  and  gamma  transport 
calculations  is  one  in  which  a  simulated  particle  is  followed  through  all  the  proc¬ 
esses  with  associated  probability  that  a  real  particle  would  follow.  Many  such 
particles  are  tracked,  and  score  their  contribution  to  any  physical  or  mathemati¬ 
cal  quantity  that  is  desired.  The  tracking  of  a  primary  particle  from  its  birth 
through  its  descendants  to  their  disappearance  from  the  given  system  (by  absorp¬ 
tion,  degradation,  or  leakage)  is  defined  as  a  history. 

A  particle  is  characterized  at  any  instant  by  its  space  position  P,  its  direction 
and  its  energy  E  (neglecting  effects  like  spin  and  interference  phenomena). 

The  first  step  is  to  determine  by  sampling  P,  and  E  of  the  virgin  particle  from 
a  normalized  source  distribution  S(P,  E,  Q). 

The  volume  anisotropic  source  generator  (VANGEN)  is  used  to  select  from  a  com¬ 
bination  of  several  types  of  spatial,  angular,  and  energy  distributions  of  the  pri¬ 
mary  particle". 

Now,  if  either  inelastic  scatterings  or  absorption  occur  in  the  course  of  tracking 
neutrons,  the  primary  Monte  Carlo  program  (ADONIS)  will  record  on  an  interac¬ 
tion  tape  the  identity  of  the  source  particle  which  is  being  tracked,  the  location 
and  the  incident  energy  of  the  neutron  at  the  interaction,  the  type  of  interaction. 
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and  the  element  with  which  the  neutron  collided.  This  is  done  for  possible  second¬ 
ary  source  generation. 

The  secondary  gamma  source  particle  generator  (GASP)  uses  the  interaction  tape 
along  with  input  data  which  characterize  the  number  and  energies  of  gammas  born 
following  an  inelastic  scattering  or  an  absorption  for  each  element  to  create  a 
secondary  source  tape  which  can  be  used  for  subsequent  Monte  Carlo  programs. 

It  should  be  noted  that  GASP  can  be  used  tc  generate  other  types  of  secondary 
source  particles  (a,  /3,  fission  products,  ...)  by  judicious  choice  of  input  data.  In 
particular,  GASP  has  been  used  to  create  a  source  for  neutrons  produced  by  fis¬ 
sion. 

An  important  feature  in  the  UNC  Monte  Carlo  programs  is  that  the  generation  of 
source  particles  has  been  arranged  in  such  a  manner  that  anyone  with  a  fundamen 
tal  knowledge  of  FORTPAN  programming  may  create  a  source  particle  generator 
complying  with  his  own  needs. 

5.1  VANGEN 

VANGEN  selects  from  a  given  combination  of  spatial,  angular,  and  energy  distri¬ 
butions  a  specified  number  of  primary  neutron  or  gamma  source  particles  and 
records  on  a  source  tape  a  number  which  identifies  the  particle,  its  position  in 
Cartesian  coordinates,  the  three  direction  cosines  of  its  initial  direction,  its  en¬ 
ergy,  and  the  region  number  in  which  the  particle  was  born. 

There  are  numerous  types  of  distributions  that  can  be  thought  of,  and  if  the  need 
arises,  a  technique  can  be  devised  to  select  from  any  of  them. 

In  the  present  programs,  we  have  restricted  ourselves  to  the  more  commonly 
used  distributions. 


5. 1. 1  Spatial  Distributions 
Volume  Source 

It  is  assumed  that  the  particles  are  uniformly  distributed  in  a  given  parallelepiped 
that  must  be  fully  cwitained  in  one  region  of  the  specified  geometrical  system.  As 
in  EZGEOM  (see  Section  6),  the  “ordinary  rectangular  parallelepiped"  must  be 
described  by  its  six  bounding  planes:  X/y,  x^,  y^,  y^,  Za, 

If  ^2)  three  random  numbers  between  0  and  1,  the  position  of  the  source 

particle  is  given  by: 


X  =  Xq;  +  (Xj3-  Xa) 

y  =  Ya  +  ^2  (y^-ya)  (42) 

Z  =  Zqj  +  ^3  (z^  ~  Zq;) 

Plane  Source 

If,  in  the  volume  source  description,  we  make  x^,  =  x^  (or  y^;^  =  y^,  or  Zq,  =  z^,  one 
obtains  the  uniform  distribution  in  a  plane  perpendicular  to  one  of  the  coordinate 
axes.  Note  that  if  a  plane  source  is  I'.sed,  it  must  be  inside  the  given  source  re¬ 
gion  and  cannot  be  coincident  with  one  of  the  six  bounding  planes  of  this  region. 

Point  Source 

If  we  now  make  xa  =  x^,  yQ,  =  Zq,  =  z^  in  the  volume  source  description  and  if 
Xq,  y^,  z^  are  the  coordinates  of  the  point  source,  then  one  gets  the  desired 
distribution. 

5.1.2  Energy  Distribution 
Method  by  Equal  Probability  Steps 

If  the  number  of  source  particles  in  energy  range  (E,  E  +  dE)  is  given  by  the  func¬ 
tion  N(E)dE,  the  input  to  VANGEN  must  be  prepared  as  follows. 
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Compute  the  probability  f{E)  by  any  available  method  of  integration: 


f(E)  = 


N(E)dE 


max 


N(E)dE 


(43) 


where  Eo  and  Eq^j^  are  the  lower  and  upper  energies  of  the  distribution  respec¬ 
tively. 

Hence,  f(Eo)  =  0,  f(Ejjjj^jj)  =  1. 

Determine  the  E^*s  such  that 

(E^) for  77  =  0,  1,  ...,  M  (44) 


where  M  is  the  number  of  equal  probability  steps  that  is  selected  (M  <100). 

The  input  will  consist  of  the  number  of  equal  probability  steps  plus  1  (M+1),  the 
width  of  the  probability  interval  (1/M)  and  the  upper  and  lower  energy  limits  of 
each  step  (i.e.,  for  step  rj,  E^_j  and  Etj). 

The  program  will  select  a  random  number,  between  0  and  1  and  take  the  inte¬ 
gral  portion  of  1  +  [^/(1/M)]  =  1  +  M^.  This  will  determine  the  energy  oin  in 
which  the  energy  of  the  source  neutron  must  be  selected.  For  example,  if  the  in¬ 
tegral  portion  of  1  +  M?  is  equal  to  ij,  the  E  $  E  ^  E  .  In  the  energy  bin  itself, 

77-1  V 

the  energy  E  is  selected  assuming  a  uniform  distribution. 

E  =E^_i  +  a(E;,-E^t)  (45) 


where 
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a 


|-[M4!  i 

jf 

M 


with  [M^]  =  integral  portion  of  M|. 

Fission  Spectrum 

The  number  of  neutrons  in  the  energy  range  (E,  E  +  dE)  is  given  by 

NiE)dE  =  A  e"“^^  sinh  dE  (46) 

where  and  A  are  constants  and  E^  <  E  <  Ey  ^  0  and  Ey  are  the  lower 

and  the  upper  limits  of  the  fission  spectrum  respectively. 

The  method  described  in  the  foregoing  section  may  be  applied.  However,  the 
program  uses  an  elegant  rejection  technique  invented  by  M.H.  Kalos®  in  order  to 
pick  E  from  the  probability  distribution  given  by  Eq.  46. 

Briefly,  an  energy  E  is  chosen  as  foEovs. 


Suppose 


Pi(x)dx= 

Xi  <  X  <  X2 

(47) 

P2(yMy  =  e"y  dy 

0  <  y  <  00 

(48) 

If 

-ax  +  (y  -bx-  c)^  <  0, 


(49) 


choose 

E  =  yx 


(50) 
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otherwise,  reject  and  repeat  from  Eq.  47. 


To  obtain  the  right  fission  distribution  and  insure  a  maximum  efficiency  for  the 
rejection  technique,  the  parameters  a,  b,  c,  y  have  been  evaluated  using  the  latest 
experimental  fit  for  the  fission  spectrum  of  prompt  neutrons. 

a  =  4.572573 
b  =  1.069179 
c  =  1.069179 
y  =  1.99675 


and  in  Eq.  46, 


Xi  = 


y 


Monoenergetic  Source 

If  the  upper  and  lower  limits  of  the  fission  spectrum  are  made  equal  to  the  desired 
energy,  one  obtains  the  monoenergetic  distribution. 

5.1.3  Angular  Distribution 
Anisotropic  Angular  Distribution 

If  the  angular  distribution  of  the  source  particles  is  symmetric  about  an  axis,no> 
one  can  divide  the  total  solid  angle  into  M  intervals,  each  interval  being  defined  by 
the  scalar  product  of  the  angle  of  the  cones  with  flo,  i.e.,  interval  t  is  defined  by 
and  Wj  such  that 


Wi_j  =  •  Oq 

Wj  =  •  flo 


(51) 
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In  each  of  these  intervals,  the  number  of  source  particles  is  given  with  probability 
Pi;  such  that 

M 

S  Pi  = 

i=l 

The  program  will  select  a  random  number  ^  between  0  and  1,  and  compare  it  to 

E  Pi> 

i=l 

where  j  <  M. 

If 

j-1  j 

E  Pi  ^  ^  ^  E  Pi» 

i=l  i=l 

then  the  direction  of  the  source  particle  is  selected  from  the  interval. 

The  direction  cosine  in  the  interval  is  then  chosen  with  equal  probability.  Thus 
W  =  Wi_i  +  |MWi-Wi_,).  (53) 

where  is  another  random  number. 

The  azimuthal  ai^le  is  also  chosen  at  random  between  0  and  2n. 

Isotropic  Angular  Distribution 

The  isotropic  case  is  obtained  by  making  M  =  l,  Wo  =  -l,Wj=+l,  and  Pj  =  1  in  the 
anisotropic  distribution. 
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^ionodi^ectional  Distribution 


In  this  case,  it  is  sufficient  to  take  M  =  1,  Wq  =  Wj  =  ±1  and  pj  =  1  in  the  anisotropic 
distribution.  Of  course,  directions  parallel  to  each  coordinate  axis  in  both  positive 
and  negative  directions  can  be  obtained  by  choosing  the  axis  direction,  fio* 

The  numbers  1,  2,  or  3  submitted  as  input  indicate  that  the  source  will  be  gener¬ 
ated  about  the  x,  y,  or  z  axis  respectively. 

5.2  GASP 

GASP  selects  the  energy  of  secondary  gammas  produced  by  neutron  inelastic  scat¬ 
tering  and  absorption  which  includes  capture,  fission,  and  production  following  a 
charged  particle  reaction.  The  position  of  these  secondary  source  gammas  is 
read  from  the  interaction  tape  and  their  direction  can  be  selected  as  in  VANGEN. 
The  output  is  written  on  the  source  tape. 

For  each  element  of  interest  (specified  by  the  integral  part  of  its  mass  number) 
for  both  absorption  and  inelastic  scattering  and  for  a  given  number  of  energy 
intervals  for  the  incident  neutrons,  the  GASP  data  input  consists  of 

1.  A  table  which  contains  the  gamma  energies  in  increasing  order, 

Eyjp  f  =  1,  2,  ...,  L  <  10. 

2.  A  table  which  lists  the  lower  and  upper  energies  of  the  neutron  bin  inter¬ 
vals,  Ejf,  k  =  0,  1,  ...,  K  <  30,  such  that  the  energy  bin  is  given  by  the 
interval  (Ej^_j,  Ej^). 

3.  A  list  of  the  expected  number  of  gammas  of  energy  Eyj^  born  in  an  inter¬ 
action  with  a  nexrtron  whose  initial  energy  lies  in  the  interval  Ej^), 

The  GASP  program  proceeds  to  select  L  random  numbers,  |j,  ^2>  •••>  •••>  II* 

The  integral  part  of  plus  a  Russian  roulette  method  for  the  fractional  por¬ 
tion  is  used  to  determine  the  number  of  source  gammas  at  energy  E^^.  For  in- 
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stance,  if  ^  =  3.5  and  is  the  corresponding  random  number,  is  compared 
to  0.5.  If  >  0.5,  then  four  gammas  are  taken  with  energy  Eyjf.  If  <  0*5,  three 
gammas  are  taken  with  energy  E^jg. 

Further,  if,  in  the  ADONIS  problem,  weights  Wjj  and  Wy  are  attached  to  the  region 
in  which  the  interaction  takes  place,  the  number  of  gammas  obtained  is  multiplied 
by  the  ratio  W^/Wy,  and  the  number  of  source  gammas  is  obtained  again  by  taking 
the  integral  fraction  of  the  product  plus  the  result  of  the  Russian  roulette  played 
with  the  fractional  portion  of  the  product. 

Gamma  production  by  radiative  neutron  capture,  fission,  and  charged  particle  re¬ 
action,  i.e.,  (N,P),  must  be  included  in  the  input  for  the  gamma  production  by  ab¬ 
sorption. 


Thus,  the  input  must  be  the  expected  number  of  gammas  born  by  neutron  absorp¬ 
tion  given  by 


(54) 


where  a  denotes  a  cross  section  and  the  superscripts  c,  F,  A,  and  P  denote  cap¬ 
ture,  fission,  absorption,  and  charged  particle  reaction,  respectively. 


5.3  VANGEN  INPUT  DESCRIPTION 

The  VANGEN  input  is  described  as  follows  (columns  73  to  80  are  for  identifica¬ 
tion): 


5.3.1  Card  1  (contains  parameter  input  for  VANGEN) 

Column  Item  Format 

I 


1-10 

15 


The  number  of  source  particles  to  be  generated. 
The  normal  direction  of  the  source:  1  =x-direc- 
tion,  2  =  y-direction,  3  =  z -direction 
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Column 

Item 

Format 

18-20 

The  number  of  the  region  where  the  soxu'ce  par¬ 
ticle  will  be  born. 

I 

23-25 

If  the  energy  distribution  is  chosen  at  equal  proba¬ 
bility  intervals,  this  entry  contains  the  number  of 
energy  intervals  +  1.  If  the  fission  spectrum  op¬ 
tion  or  monoenergetic  option  is  desired,  the  entry 
is  0. 

I 

26-37 

The  entry  for  the  upper  energy  limit  of  the  fission 
spectrum  desired,  Eu  (0  if  this  option  is  not  de¬ 
sired). 

E 

38-49 

Note: 

The  entry  for  the  lower  energy  limit  of  the  fission 
spectrum  desired,  El  (0  if  this  option  is  not  de¬ 
sired). 

E 

For  monoenergctic  problems,  and  entries  in  the  columns  23 

through  ?5  should  be  0. 


5.3.2  Card  2  (contains  the  desired  spatial  distribution) 


Colunrn 


Item 


Format 


1-12 

12-24 

25>.36 

37-48 

49-60 

61-72 

Notes: 


The  X(y  bounding  plane  of  the  source  volume. 
The  bounding  plane  of  the  source  volume. 
The  yq  bounding  plane  of  the  source  volume. 
The  bounding  plane  of  the  source  volume. 
The  Zq  bounding  plane  of  the  source  volume. 
The  z^  bounding  plane  of  the  source  volume. 


6E12 


1.  If  Xq  =  x^,  or  ya  =  y^t_2E  ^ot  =  plane  source  distribution  is  ob¬ 
tained. 

2.  If  Xfl  =  x^,  yo,  =  y^  and  Zq  =  z^g,  point  source  distribution  is  obtained. 


At  the  present  time,  the  spatial  distribution  must  be  contained  in  one  and  only  one 
region. 


The  next  set  of  cards  contains  information  on  the  energy  distribution  in  equal 
probability  steps  if  desired.  (Use  only  if  columns  23  to  25  on  card  1  are  nonzero.) 
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The  first  entry  is  the  probability  step.  It  is  then  followed  by  the  energies  (in  ev) 
in  increasing  order,  (These  energies  are  the  limits  of  energy  bins.)  The  number 
of  energies  is  one  more  than  the  number  of  bins.  There  are  six  entries  per  card. 
The  format  is  6E12.5. 

Three  types  of  cards  describe  the  angular  distribution: 

1.  The  Bin  Information  Card 

Contains  the  number  of  angular  bins  +  1  in  columns  3-4. 

2.  The  Angular  Bin  Card 

These  cards  contain  the  cosine  limits  of  each  angular  bin.  There  are 
seven  entries  per  card  and  the  format  is  7E10.4.  [The  lower  limit  of 
the  (i+1)^^  bin  is  the  upper  limit  of  bin.] 

3.  The  Probability  Card 

The  first  entry  is  0  and  each  entry  contains  the  corresponding  probability 
for  each  angular  bin  starting  with  the  first  angular  bin.  There  are  seven 
entries  per  card  and  the  format  is  7E10,4. 

Note: 

If  an  isotropic  source  in  the  forward  half  plane  is  desired,  only  one  angu¬ 
lar  bin  is  necessary.  The  input  then  is: 

Card  type  1  2  Two  cosines  define  bin  1 

Card  type  2  0.0  1.0  The  lower  and  upper  limits  of  the  bin 

Card  type  3  0  1.0  The  probability  that  a  particle  will  be  born 

in  bin  1  is  1. 
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5.4  GASP  INPUT  DESCRIPTION 


The  GASP  input  is  described  as  follows: 

5.4.1  Element  Parameter  Card 

This  card  contains  the  number  of  elements  used  in  the  associated  ADONIS  problem. 
The  format  is  15. 

5.4.2  Element  Input  Data  (for  each  element) 

1.  Element  Parameter  Card 


Columa 

Item 

Format 

-•  r 

Integral  part  of  the  atomic  mass  (identification  of 
element  same  as  that  specified  in  DATORG  input) 

6-10 

La  =  number  of  absorption  gamma  output  energies 
(«:10)  (0  if  no  capture  y  produced) 

I 

11-15 

=  number  of  corresponding  incident  neutron  en¬ 
ergies  (<30).  This  equals  the  number  of  neu¬ 
tron  energfy  bins  +  1  (0  for  no  y  capture). 

I 

16-20 

Lj  =  number  of  inelastic  y  out  energies  (<10) 

I 

21-25 

Ki  =  number  of  corresponding  incident  neutron 
energies  (^30).  This  equals  the  number  of 
neutron  energy  bins  +  1. 

I 

2.  Capture  y  Energy  Cards  (as  many  as  necessary) 

The  energies  are  given  in  ascending  order  in  ev.  The  number  of  entries 
per  card  is  six  and  the  format  is  6E12.6.  The  total  number  of  entries  is 
equal  to  L^. 

3.  Capture  Neutron  Energy  Cards  (as  many  as  necessary) 

Every  two  entries  represent  a  neutron  energy  bin.  The  entries  are  in 
ascending  order  and  in  ev.  There  are  six  entries  per  card  and  the  format 
is  6E12.6. 
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4.  Capture  Number  Cards 


Each  entry  represents  the  number  of  gammas  produced  at  each  of  the 
specified  y  energies.  The'^e  numbers  are  entered  for  each  neutron 
energy  bin.  The  first  entries  are  0  reprv^senting  the  0^^  bin. 

Start  a  new  card  for  each  neutron  bin.  There  are  six  entries  per  card 
and  the  format  is  6E12.6. 

Repeat  items  2,  3,  and  4  for  inelastic  scattering  input,  if  any. 

5.4.3  Weights  for  Neutron  and  Gamma  Problem  Input  (See  Section  7.3.1) 

1.  Number  of  Regions  Card 

Enter  the  number  of  regions  used  in  the  associated  neutron  and  gamma 
problems  (15  format). 

2.  Enter  the  weight  for  each  region  used  in  the  neutron  problem  (for¬ 
mat  5E14.8).  (Note;  there  is  only  one  entry  per  region.  At  the  present 
time  one  cannot  use  energy  importance  in  the  neutron  problem  and  run 
the  secondary  y-ADONIS.  It  is  anticipated  that  this  restriction  will  be 
removed  in  the  near  future.) 

3.  Enter  the  weight  for  each  region  to  be  used  in  the  secondary  gamma 
problem  (format  5E14.8). 

5.4.4  Angular  Distribution 

The  angular  distribution  input  is  the  same  as  for  VANGEN. 

5.5  DESCRIPTION  OF  THE  SOURCE  TAPE 

The  source  tape  which  is  generated  by  both  VANGEN  and  GASP  must  conform  to 
the  specifications  required  by  the  ADONIS  progi'am.  Particles  must  be  generated 
in  groups  of  100.  Each  particle  is  described  by  nine  numbers: 
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1  the  particle  number 

2-4  the  location  of  the  particle  in  3 -space  (its  x-y-z  coordinates) 

5-7  the  three  direction  cosines  the  particle  assumes  initially 

8  the  initial  energy  of  the  particle 

9  the  region  number  in  which  the  particle  is  born. 

With  the  exception  of  the  particle  number  and  the  region  number,  the  entries  are 
in  floating  point  notation,  and  are  generated  by  a  FORTRAN  routine  using  the 
^*WRITE  TAPE**  statement  for  100  particles,  with  nine  entries  for  one  particle 
preceding  the  nine  entries  for  the  next.  Finally,  to  signal  the  end  of  the  source 
particles  that  one  wants  to  generate,  one  additional  particle  whose  particle  num¬ 
ber  is  -1  is  necessary. 

5.6  SOURCE  INPUT 

The  source  input  consists  of  either  VANGEN  or  GASP  input  forms.  See  following 
pages. 
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6.  THE  ORGANIZATION  OF  THE  GEOMETRY  DATA  (EZGEOM) 


Monte  Carlo  programs  track  simulated  particles  through  a  specified  geometrical 
configuration,  undergoing  all  the  interactions  that  an  actual  particle  is  expected  to 
undergo.  The  particle  flight  may  pass  through  media  having  different  properties. 
Therefore,  it  is  of  primary  importance  to  know  at  all  times  in  which  medium  the 
particle  is,  when  the  particle  leaves  the  medium,  and  which  n.edium  it  will  be  en¬ 
tering.  Different  media  may  be  described  as  containing  given  material  composi¬ 
tions  in  three-dimensional  geometrical  figures. 

However,  it  is  necessary  to  decompose  these  figures  into  particular  geometrical 
forms  which  can  be  handled  by  the  ADONIS  program,  i.e.,  nonintersecting  rectan¬ 
gular  parallelepipeds,  right  cylinders,  spheres,  and  right  wedges. 

The  EZGEOM  program  processes  the  input  which  consists  of  the  description  of 
the  desired  configuration  by  means  of  these  elementary  geometrical  forms  and 
stores  the  information  into  tables  for  rapid  data  access  in  subsequent  Monte  Carlo 
programs. 

6.1  DEFINITIONS  AND  NOTATION 

The  three-space  coordinate  system  used  in  describing  the  geometry  is  a  Cartesian 
right-handed  coordinate  system.  The  x,  y,  z  directions  will  be  described  as  Xj, 

X2,  X3. 
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The  entire  three-dimensional  space  is  initially  decomposed  into  a  set  of  rectangu¬ 
lar  parallelepipeds  whose  bounding  planes  are  parallel  to  the  coordinate  axes. 
These  parallelepipeds  are  caller  ordinary  regions  and  are  described  by  the  co¬ 
ordinates  of  their  six  bounding  planes:  xp,  x^  with  i  =  1,  2,  3  and  x^  <  x^. 

Any  elementary  geometrical  form,  which  is  completely  described  in  the  input  by  a 
vertex  V  (x^,  X2,  X3)  and  a  set  of  vectors  and  scalars  is  called  a  nonordinary  re¬ 
gion.  Nonordinary  regions  presently  handled  by  EZGEOM  and  geometry  tracking 
routines  are  spheres,  right  cylinders,  rectangular  parallelepipeds,  and  right 
wedges.  In  the  near  future,  right  frusta  of  cones  will  be  liandled  by  these  routines. 

Due  to  storage  restrictions  and  word  length  requirements,  the  following  rules 
must  be  observed  when  preparing  E liGEOM  input. 

1.  Nonordinary  regions  must  be  wholly  contained  in  either  an  ordinary  or  a 
nonordinary  region.  Due  to  possible  round-off  errors,  it  is  suggested 
that  no  surface  of  a  contained  region  be  closer  to  any  surface  of  the  con¬ 
taining  region  than  10"®  times  the  maximum  coordinate  found  on  the  sur¬ 
face  of  the  containing  region.  Note  that  the  geometry  is  first  decomposed 
into  ordinary  regions.  Hence,  a  nonordinary  region  may  contain  a  set  of 
nonordinary  regions,  but  must  eventually  be  itself  part  of  a  set  of  regions 
which  are  wholly  contained  in  an  ordinary  region. 

2.  The  total  number  of  ordinary  and  nonordinary  regions  must  be  less  than 
or  equal  to  200. 

3.  At  most,  127  nonordinary  regions  may  be  used. 

4.  A  nonordinary  region  I  must  be  wholly  contained  in  another  region  E.  Let 
region  E  be  defined  as  the  region  external  to  I.  At  the  most,  seven  such 
regions  may  have  E  as  an  external  region.  We  will  say  that  any  region 
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may  see  only  seven  regions.  In  turn,  each  of  these  regions  may  see,  at 
most,  seven  regions. 

5.  The  geometry  input  must  be  submitted  as  follows: 

First,  define  the  ordinary  regions  and  number  them  from  1  to  M,  M  being 
the  number  of  ordinary  regions. 

Then,  define  the  nonordinary  regions  (by  geometrical  category  and  co¬ 
ordinate  data)  and  number  them  starting  from  M  +  1. 

6.  The  data  required  to  describe  the  nonordinary  regions  for  the  input  are: 

a.  Spheres 

Specify  the  vertex  V  at  the  center  and  the  scalar  R  denoting  the  radius. 

b.  Right  Cylinders 

Specify  the  vertex  V  at  the  center  of  one  base,  a  height  vector,  H,  ex¬ 
pressed  in  terms  of  its  components,  and  a  scalar  R  denoting  the  base 
radius. 

c.  Rectangular  Parallelepipeds 

Specify  the  vertex  V  at  one  of  the  corners  and  a  set  of  three  mutually 
perpendicular  vectors,  aj,  representing  the  height,  width,  and  length 
of  the  parallelepiped  respectively.  Note  that  the  bounding  planes  need 
not  be  parallel  to  the  coordinate  axes. 

d.  Right  Wedges 

The  input  is  the  same  as  for  the  nonordinary  rectangular  parallele¬ 
pipeds.  However,  the  first  two  vectors  describe  the  two  legs  of  the 
right  triangle  of  the  wedge. 
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e.  Right  Frusta  of  Cones 


Specify  a  vertex  V  at  the  center  of  the  lower  base,  the  height  vector,  H, 
expressed  in  terms  of  its  components,  and  two  scalars,  R,  and  R,,  de¬ 
noting  the  radii  of  the  lower  and  upper  bases  respectively.  (The  frusta 
of  cones  is  e.xpected  to  be  included  in  the  geometry  description  in  the 
near  future.) 

6.2  TECHNICAL  DESCRIPTION 

The  EZGEOM  routine  considers  first  the  ordinary  regions  and  proceeds  to  find 
the  regions  adjacent  to  each  side  of  the  ordinary  parallelepipeds,  starting  with 
region  number  1. 

cy  3 

^i  m'  ^  m  coordinates  of  the  six  bounding  planes  of  the  ordinary  re¬ 

gion,  m  (1  <  m  M,  M  being  the  total  number  of  ordinary  regions).  The  routine 

ct  Q 

starts  with  Xj^j  and  compares  it  to  xj  for  m  =  2,  3,  ...,  M. 

Suppose  that  xJ^j  =  ^,mj  then  the  xf  face  of  region  1  and  the  face  of  region  m 
are  on  the  same  plane.  This  case  is  illustrated  by  Fig.  4, 

It  is  necessary  to  determine  whether  or  not  the  two  faces  have  an  area  in  common. 
Hence,  the  program  proceeds  to  test: 

1.  x^j:  xf^m-  If  greater  than  or  equal,  there  is  no  overlap  and  the  routine 
considers  the  next  region  (Fig.  4  and  case  a  of  Fig.  5).  If  less  than,  then 
the  second  test  is  tried. 

2.  j:  If  less  than  or  equal,  there  is  no  overlap  and  the  routine  con¬ 

siders  the  next  region  (case  a.  Fig.  5).  If  greater  than,  then  the  third 
test  is  tried. 
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3.  x^jt  X3  m.  If  greater  than  or  equal,  there  is  no  overlap  (case  b,  Fig.  5) 
and  the  comparison  to  the  m  +  1  region  is  considered,  if  less  than,  then 
the  fourth  and  last  test  is  tried. 

4.  x£i;  If  less  than  or  equal,  there  is  no  overlap  (case  b.  Fig.  5)  and 

the  comparison  to  the  m  +  1  region  is  considered.  If  greater  tiian,  then 
there  is  an  area  in  common  (case  c.  Fig.  5).  In  this  case,  the  region  m 

is  stored  in  a  table  in  which  all  regions  adjacent  to  region  1  are  recorded. 

Once  the  comparison  of  x^j  to  x^^m  ^  Is  ccrrpleted,  the  number  of  regions 

contained  in  the  table  is  tested. 

1.  If  there  is  no  adjacent  region,  the  coordinate  x^j  is  an  extremity  for  the 
geometry  and  le  region  adjacent  to  region  1  on  the  x®  face  is  itself. 

This  coordinate  is  saved.  Then,  if  in  the  course  of  the  investigation  of 
the  X®  coordinate  of  other  regions  another  extremity  is  found  which  is 
not  equal  to  x^j,  it  means  that  the  geometry  has  not  been  properly  de¬ 
fined.  This  diagnostic  is  printed,  and  the  program  is  stopped. 

2.  If  there  is  only  one  region,  the  region  number  is  placed  in  the  proper 
table  as  the  region  adjacent  to  region  1  on  the  x®  face,  f'hirther,  a  flag 
is  set  to  denote  that  region  1  is  adjacent  to  that  region,  say  region  m. 
Thus,  when  the  x^  face  of  region  m  will  be  investigated,  there  is  no  need 
to  search  for  regions  in  contact  with  it  if  their  region  number  is  less 
than  m.  The  procedure  is  followed  for  regions  whose  number  is  greater 
than  m  only. 

ot 

3.  If  there  is  more  than  one  region  in  contact  with  the  Xj^j^  face  of  a  given 

region,  the  face  is  defined  as  a  complex  surface.  The  EZGEOM 
program  proceeds  to  decompose  the  complex  surface  into  a  grid  of 
elementary  rectangles  and  assign  to  each  of  them  the  corresponding  re¬ 
gion  number. 
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Let  us  consider  the  complex  surface  oriented  such  that  the  normal  to  that  plane  is 
in  the  positive  direction.  The  horizontal  and  the  vertical  directions  are  deter¬ 
mined,  using  the  right-handed  method,  i.e,,  if  the  X2-axis  is  normal  to  the  complex 
surface  a.'^  in  the  case  shown  in  Fig.  6,  then  the  Xj  and  X3  axes  are  the  horizontal 
and  vertical  directions  respectively. 

The  first  step  performed  by  the  routine  is  to  extract  the  ^  coordinates  in  the  hori¬ 
zontal  direction  of  all  adjacent  regions  and  place  them  in  increasing  order  remov¬ 
ing  duplicated  coordinates,  say  Hj,  H2,  ...,  Hj^ax* 

The  second  step  is  to  test  if  greater  than  or  equal  to  xf  of  the  region  in¬ 
vestigated.  If  it  is,  xf  replaces  Hj^ax  table.  If  it  is  less,  this  means  that 

there  is  an  area  on  the  right  side  of  the  face  which  has  not  been  accounted  for.  A 
diagnostic  is  printed  pointing  out  that  a  hole  exists  and  the  program  is  stopoed. 

The  same  sequence  of  operations  is  performed  for  the  j3  coordinates  in  the  vertical 
direction,  say  Vj,  V2,  ...,  V^ax* 

These  horizontal  and  vertical  sequences  determine  a  two-dimensional  matrix.  The 
complex  surface  is  thus  decomposed  in  a  grid  of  rectangles.  The  EZGEOM  pro¬ 
gram  proceeds  to  place  in  each  of  these  rectangles  the  corresponding  region  num¬ 
ber.  It  must  be  noted  that  several  rectangles  may  belong  to  a  given  region.  This 
is  done  as  follows:  let  x^m,  be  the  coordinates  in  the  horizontal  direction  of 
the  adjacent  region  m.  The  program  determines  the  num.ber  of  intervals  Hi  con¬ 
tained  in  [x®,x^].  Clearly,  there  must  be  horizontal  coordinates  in  the  horizontal 
sequence  which  equal  to  xf  and  x^.  If  not,  either  a  hole  or  an  overlap  exists  and  a 
diagnostic  is  printed  giving  the  region  number  tnder  consideration,  the  plane  being 
searched,  the  table  of  adjacent  regions  found,  and  the  horizontal  and  vertical 
tables.  From  this  iriformation,  it  is  possible  to  construct  the  complex  surface  and 
find  the  error. 
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Illustration  of  tlie  EZGEOM  decomposition  of  a  complex 


If  no  error  has  occurred,  the  data  gathered  for  the  complex  surface  are  placed  in 
Table  3  which  is  described  in  Section  6.3. 

The  same  procedure  is  repeated  for  each  face  of  each  region. 

When  all  ordinary  regions  have  been  processed,  the  input  cards  for  the  nonordi¬ 
nary  regions  are  read  and  stored  according  to  the  format  described  in  Section  6.3. 

At  the  present  time,  there  is  no  routine  built  into  EZGEOM  which  verifies  that  the 
submitted  input  for  the  nonordinary  region  is  properly  defined.  Consequently,  ex¬ 
treme  care  should  be  taken  in  preparing  input  data  for  nonordinary  regions. 

In  conclusion,  it  is  believed  that  EZGEOM  will  find  most,  if  not  all,  errors  made 
in  the  input  data  for  ordinary  geometry.  It  is  possible  for  EZGEOM  to  fail  to  give 
a  diagnostic  when  the  prc^rammed  dimensions  of  a  table  for  complex  surfaces  are 
exceeded.  This  may  occur  when  a  complex  surface  is  defined  by  a  grid  of  more 
than  20  x  20  rectangles  and  the  investigated  region  is  located  near  the  end  of  the 
ordinary  region  listing.  This  is  an  instance  when  an  incorrect  edit  can  be  ob¬ 
tained.  In  this  case,  all  the  regions  adjacent  to  region  1  will  oe  called  region  1, 
showing  that  the  geometry  was  handled  incorrectly.  In  general,  the  situation  is 
improved  by  assigning  a  very  low  region  number  to  that  region. 

6.3  DESCRIPTION  OF  THE  INPUT-OUTPUT 

6.3.1  Input 

1.  The  first  card  contains  the  number  of  ordinary  regions,  M,  and  the  num¬ 
ber  of  nonordinary  regions,  N  (N<127).  Note  that  M  +  N  <  200. 

2.  The  next  M  cards  contain  the  coordinates  of  the  six  bounding  planes  of 
each  region:  x®,  i  =  1,  2,  3;  x^*  <  x^.  Note  that  the  cards  are  placed 
in  the  same  order  ae  the  region  numbers. 
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3.  The  input  for  each  of  the  N  nonordinary  regions  has  been  described  in  the 
definitions  (Section  6.1).  The  order  in  which  the  cards  are  placed  is  also 
according  to  the  region  numbers  (M  +  1  to  N  +  M). 

6.3.2  Output 

The  output  edit  gives  a  complete  description  of  each  region  and  is  self-explanatory. 

EZGEOM  generates  three  tables  which  completely  describe  the  geometry  specifi¬ 
cations  for  the  problem.  The  data  in  these  tables  are  arranged  so  that  they  are 
rapdily  accessible  to  the  Monte  Carlo  tracldng  routines. 

Table  1 

There  is  an  entry  in  Table  1  corresponding  to  each  region  of  the  problem.  (The 
entry  corresponds  to  the  i^^  region.)  Each  entry  to  Table  1  is  packed  into 
35  bits  as  follows. 

1.  The  first  three  high  order  bits  contain  the  number  of  internal  regions 
that  the  region  sees  (see  Section  6.1).  The  maximum  number  of  such 
regions  for  any  given  region  is  se  v^en. 

2.  The  next  three  bits  describe  the  i^pe  of  region 
0  =  ordinary  parallelepiped 

1  =  sphere 

2  =  cylinder 

3  =  wedge 

4  =  nonordinary  parallelepiped 

5  =  frusta  of  cone. 

3.  The  next  eight  bits  contain  the  number  of  the  region  external  to  the  re¬ 
gion.  This  only  applies  to  nonordinary  regions,  types  1-5). 

4.  The  next  seven  bits  contain  a  location  relative  to  Table  2  where  the  re¬ 
gions  whose  external  region  is  the  region  considered  are  listed.  The 
entry  is  "’.ero  when  there  are  no  such  regions. 
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5.  The  last  14  bits  contain  the  location  relative  to  Table  3  where  the  region 


data  are  listed. 


Table  2 

This  table  contains  an  entry  for  each  of  the  nonordinary  regions  (types  1-5),  and 
these  entries  are  grouped  corresponding  to  their  external  region.  Each  entry  is 
packed  into  25  bits  as  follows. 

1.  The  first  three  bits  contain  the  type  of  the  region  (same  as  Table  1,  item  2 
but  0  type  is  excluded). 

2.  The  next  eight  bits  contain  the  region  number  of  the  nonordinary  region. 

3.  The  next  14  bits  contain  the  location  relative  to  Table  3  where  the  region 
data  are  listed. 

Table  3 

This  table  contains  the  parameters  that  define  each  region.  The  first  location  is 
referenced  from  Table  1  or  Table  2. 

1.  Ordinary  Parallelepipeds 

The  first  six  entries  contain  the  six  bounding  planes  (x®,  xf,  i  =  1,  2,  3). 
The  next  six  entries  contain  the  adjacent  region  numbers  for  each  of  the 
corresponding  six  sides. 

a.  If  a  side  sees  more  than  one  region  this  side  is  defined  as  a  complex 
surface.  The  entry  contains  the  negative  of  the  location  relative  to 
Table  3  where  the  data  for  the  complex  surface  may  be  found.  At  this 
location  relative  to  Table  3  the  following  information  is  recorded: 

(1)  A  word  containing  the  number  of  horizontal  and  the  number  of  ver¬ 
tical  intervals  for  the  complex  surface  packed  18  bits  apart. 

(2)  Following  this  word  appear  the  upper  limit  of  each  of  the  horizon¬ 
tal  intervals. 
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(3)  Then  the  upper  limits  of  each  of  the  vertical  intervals  appear. 

(4)  Finally,  packed  four  to  a  word,  nine  bits  apart,  are  the  regions  on 
the  complex  surface.  The  region  numbers  along  the  horizontal, 
from  the  first  to  the  last  horizontal  interval  are  listed  for  the 
first  vertical  interval;  then  for  the  second  vertical  interval,  etc. 
as  shown  in  Fig.  6. 

b.  If  the  side  sees  only  one  region  then  the  region  number  of  this  region 
appears. 

2.  Nonordinary  Regions  (Types  1-5) 

a.  Spheres 

The  first  three  entries  contain  the  coordinates  of  the  center  of  the 
sphere.  The  next  entry  is  the  square  of  the  sphere  radius. 

b.  Cylinders 

The  first  three  entries  contain  the  coordinates  of  the  center  of  the 
base.  The  next  entry  contains  the  square  of  the  base  radius.  The 
next  three  entries  contain  the  x-y-'^  components  of  the  height  vector 
of  the  cylinder. 

c.  Nonordinary  Par>\llelepipeds 

The  first  three  entr'es  contain  the  coordinates  of  the  vertex  of  the 
parallelepiped.  The  next  nine  entries  contain  the  x-y-z  components 
of  the  three  vectors  which  define  the  parallelepiped, 

c.  Wedges 

The  description  for  the  wedge  is  the  same  as  the  nonordinary  parallele¬ 
piped,  the  two  first  vectors  defining  the  legs  of  the  wedge  right  tri¬ 
angle. 

e.  Right  Frusta  of  Cones 

The  first  three  entries  contain  the  coordinates  of  the  center  of  the  low¬ 
er  base.  The  next  two  entries  contain  fhe  square  of  the  radii  of  the 
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lower  and  upper  bases  respectively.  The  following  three  entries  con¬ 
tain  the  x-y-z  components  of  the  height  vector. 

6.3,3  Input  Sheets  for  EZGEOM 

Input  sheets  for  EZGEOM  are  shown  on  the  following  pages. 
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BOUNDING  PLANE  CARDS  FOR  ORDINARY  RECTANGULAR  PARALLELEPIPEDS  (ONE  CARD  PER  REGION) 


'This  number  should  be  the  number  of  the  pertinent  region. 


EZGEOM  INPUT  {CONTINUED) 


EXGEOM  INPUT  (CONTINUED) 


EZGEOM  INPUT  (CONTINUED) 


Repeat  above  3  carets  for  each  wedge 


EZGEOM  INPUT  (CONTINUED) 


Repe^it  cibove  2  cards  for  each  fnishim 
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7.  THE  TRACONG  AND  SCORING  MONTE  CARLO  PROGRAM  -  ADONIS 


Having  appropriate  output  tapes  of  DATORG,  VANGEN  or  GASP,  and  EZGEOM, 
we  are  ready  to  proceed  with  the  transport  game’^  for  the  solution  of  the  Boltzmann 
integral  equations.  The  transport  game  is  basically  an  analog  of  the  physical 
transport  of  a  neutron  or  a  gamma  ray.  However,  in  most  problems,  it  is  neces¬ 
sary  to  modify  the  game  to  obtain  the  desired  result  efficiently. 

At  each  event  to  be  considered  in  the  game,  a  weight  is  introduced.  The  purpose 
of  carrying  along  the  weight  is  to  enable  the  consideration  of  altered  games  of 
chance  which  lead  to  appreciable  variance  reductions.  At  the  present  time,  the 
ADONIS  program  employs  two  methods  of  variance  reduction:  importance  sam¬ 
pling,  and  statistical  estimations  in  the  special  case  of  flux  evaluation  in  small 
detector  regions. 

As  the  transport  game  is  played,  a  variety  of  information,  including  the  following, 
is  recorded. 

1.  The  average  flux  as  a  function  of  energy  group  in  each  ordinary  and  non¬ 
ordinary  region  along  with  its  standard  deviation.  If  any  quantity  pro¬ 
portional  to  the  flux  is  desired,  it  can  be  computed  using  appropriate 
response  functions. 

2.  A  neutron  interaction  tape  for  future  secondary  source  emission  process¬ 
ing  . 
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3.  A  tape  containing  information  on  transmitted  particles  for  future  process¬ 
ing. 

4.  The  number  of  absorptions  per  element  and  per  region,  births  and 
deaths  due  to  importance  sampling  in  each  region,  and  energy  degrada¬ 
tions  below  an  optional  energy  cutoff. 

5.  The  average  flux  as  a  function  of  energy  at  specified  detector  points  along 
with  its  standard  deviation  (optional,  for  neutron  problems  only). 

We  turn  our  attention  to  the  description  of  the  transport  game  played  by  ADONIS. 

7.1  THE  TRACKING  PROCESS 

The  ADONIS  program  uses  the  source  particle  tape  created  by  VANGEN,  GASP, 
or  a  similar  source  generator  program  to  start  particle  histories.  The  position, 

P,  the  energy,  E,  and  the  direction,  ,  of  the  source  particle,  and  R^,  the  region 
in  which  the  particle  is  born,  are  read  off  the  tape. 

The  first  problem  that  must  be  solved  is  to  determine  the  position  of  the  next 
collision  and  the  region  in  which  that  collision  will  occur.  It  can  be  shown  that 
the  number  of  mean  free  paths.  A,  that  a  neutron  is  to  travel  before  interacting 
must  be  selected  from  a  probability  distribution  function  e"^  (using  the  definition 
of  the  transport  kernel  in  the  Boltzmann  integral  equation).  Then,  the  distance 
from  P  along  ^  to  the  first  collision  location  P'  is  given  by 

n 

A  -  2  MT(E,Rk)  Sk  =  0  (55) 

k=i 

where  p'p(E,Rk)  =  the  total  cross  section  at  energy  E  of  the  composition  in 

region  Rj^ 

Sk  =  the  distance  along  the  particle  flight  path  in  region  Rk* 

Rf,  Rj+j,  ...,  Rk>  •••>  are  the  regions  which  are  seen  successively  along  the 
flight  path. 
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Further,  it  can  be  shown  that,  if  ^  is  a  random  number  in  the  interval  (0,  1),  the 
probability  distribution  for  -In  ^  is  e""^. 

Therefore, 

A  =  -In  ^  (56) 

The  techniques  used  by  the  ADONIS  program  to  solve  the  problem  now  can  be 
discussed. 

1.  A  random  number,  4,  is  selected  using  the  subroutine  RANDM.  This 

subroutine  is  written  in  FORTRAN  language  and  generates  equi-distribut- 
ed  pseudorandom  numbers  on  the  interval  (0,  1)  compatible  with  the  word 
length  of  the  computing  machine  used.  For  the  CDC-1604-A,  the  period 
of  the  random  number  sequence  is  Xq  =  5'^  and  the  starting  number  is 
usually  an  odd  prime  number.  The  successor,  of  a  random  number, 

is  obtained  by  taking  the  product  Xo^n>  rejecting  the  first  half  of  the 
significant  figures  of  the  product  and  floating  the  remainder. 

2.  The  number  of  mean  free  paths  A  is  computed  using  Eq.  56. 

3.  The  cross  section  M'j’(E,Ri)  of  the  composition  in  the  region  where  the 
particle  is  born  is  computed  by  linear  interpolation  of  the  cross  seriions 
generated  by  DATORG  at  the  output  energies  defined  in  the  GENOA  rou¬ 
tine.  The  subroutine  Data  Retrieval,  entry  1  (DR,1)  computes  the  inter¬ 
polation  factor,  f.  The  subroutine  Data  Retrieval,  entry  2  (DR, 2)  com¬ 
putes  M'p(E,Ri)  using  the  interpolation  factor  obtained  from  (DR,1)  and 
the  table  of  total  cross  sections  for  the  composition  in  Rj.  Note  that  in 
going  through  regions  containing  different  compositions  without  change 
in  particle  energy,  only  (DR,2)  is  called  since  the  interpolation  'actor  is 
unchanged. 

4.  Given  the  location  P  and  the  direction  £2  in  a  given  region,  R^,  the  sub¬ 
routine  Geometry  1  (Gl)  determines  how  far  it  is  to  the  next  region  and 
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what  surface  the  particle  will  cross  if  it  does  travel  that  far,  using  the 
geometry  tables  created  by  EZGEOM  (Section  6)  and  special  routines 
for  each  type  of  region,  G1  (TYP),  where  TYP  =  ORP,  SPH,  CYL,  WED, 
RPP.  These  denote  ordinary  region,  sphere,  cylinder,  wedge,  or  non- 
ordinary  parallelepiped,  respectively. 

The  number  of  regions,  mg,  wholly  contained  in  a  given  region  Ri  is 
obtained  from  Table  1, 

If  mg  =  0,  Gl(TYP)  and  R^  determine  Sj.  The  intersectio;:  P'  of  the 
flight  path  with  the  appropriate  edge  of  R^  is  computed*  and 

Si  =  I  P'-PI. 

If  m  *  0,  G1  determines  whether  or  not  tlie  flight  path  intersects  each 
B 

of  the  mg  regions  enumerated  in  Table  2  using  Gl(TYP).  If  there  is  no 
intersection  with  any  of  these  regions,  the  procedure  is  similar  to  the 
case  mg  =  0.  If  there  are  one  or  more  intersections,  the  distances  from 
P  tc  the  intersection  points  are  computed  and  the  shortest  distance  is 
then  retained. 

5.  The  next  task  is  to  determine  whether  or  not  the  particle  makes  a  col¬ 
lision  before  traveling  the  distance  S^. 

If  A  -  p.j,(E,Rp  >  0,  the  particle  does  cross  into  the  next  region.  A  subroutine 
called  Geometry  2  (G2)  is  now  used  to  determine  the  region,  Rj+j,  and  the  particle 
enters  using  the  information  gathered  in  Tables  2,  3,  and  Gl.  Then  the  remaining 
number  of  mean  free  paths  is  determined  by  replacing  A  by 


A'  =  A  -  iUT(E,Hi)  Sj  (57) 

and  the  procedure  is  repeated  by  computing  M-p(E,R^^j)  and  going  to  Gl.  In  other 


*The  method  used  in  these  computations  has  been  described  in  unpublished  memos 
by  J.  Certaine.®>®>^“ 
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words,  the  particle  is  tracked  into  the  next  region.  There  are,  however,  two  im¬ 
portant  tests  that  are  performed  before  tracking  in  the  next  region  is  started. 

First,  the  program  determines  whether  tl*e  region  in  which  the  particle  enters 
is  a  reflecting  or  transmitting  region.  * 

If  this  is  the  case,  the  tracking  of  the  particle  is  terminated  and  the  proper  reflect¬ 
ing  or  transmitting  counters  are  updated. 

If  the  region  is  neither  reflecting  nor  transmitting,  a  second  test  is  performed. 

As  pointed  out  above,  the  consideration  of  variance  reduction  techniques  may  lead 
to  the  introduction  of  weighting  functions  which  are  position  and  energy  dependent. 


At  the  present  time,  ADONIS  considers  weighting  functions  piecewise  constant  in 
each  region  and  in  given  energy  intervals,  W(E,R)  (Section  7.3,1).  The  program 
compares  the  ratio  of  weights  at  energy  E  in  region  and  R^+j. 


*It  is  necessary  to  terminate  the  particle  tracking  upon  leaving  tne  environment 
of  the  geometrical  configuration  of  interest.  Three  devices  have  been  built  into 
the  program  to  handle  this  situation. 

1.  If  the  “importance"  of  a  particle  in  a  particular  region  is  made  sufficient¬ 
ly  low,  the  particle  entering  that  region  will  be  killed  (see  Section  7.3). 

2.  If  a  fictitious  macroscopic  absorption  cross  section  sufficiently  large 
is  attached  to  a  given  region,  incident  particles  v/ill  be  absorbed. 

3.  A  very  rapid  way  to  terminate  particle  tracking  is  accomplished  by  the 
introduction  of  tw'O  special  types  of  regions:  the  reflecting  region  (re¬ 
quired  to  be  region  1)  and  the  transmitting  region  (which  can  be  any 
region)  which  have  the  following  properties: 

a.  Particles  incident  on  both  regions  are  terminated  and  counted  sepa¬ 
rately. 

b.  For  transmitted  particles  only,  the  location  at  which  the  particle 
path  flight  enters  the  region,  the  energy,  the  direction,  and  the  orig¬ 
inal  history  number  of  the  particle  and  the  region  from  which  the 
particle  has  come  are  recorded  on  a  tape  in  the  same  order  as  speci¬ 
fied  for  the  source  tape  (Section  5), 
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(58) 


W(E,Ri) 

W(E,Ri^j) 


=  N  +  F 


where  N  and  F  are  the  integral  and  fractional  portions  of  the  quotient,  respective¬ 
ly. 

A  random  number  |  is  selected  and  compared  to  F;  then,  if  I  <  F,  N  is  replaced 
by  N+1.  If  k  >  F,  N  is  left  unchanged.  N  is  by  definition  the  number  of  particles 
creat?*d.  vj  =  0,  the  tracking  is  terminated  and  the  counter  registering  the  num¬ 
ber  of  deaths  by  importance  sampling  is  increased  by  1.  If  N  =  1,  no  additional 
particle  has  been  created  and  the  tracking  of  the  particle  is  continued.  Note  that 
if  W(E,R)  =  1  for  :  .1  E  and  R,  this  means  that  the  transport  game  is  modeled 
directly  after  the  physical  problem  of  particle  transport.  If  N  >  1,  the  phase- 
space  coordinates  of  the  N-1  additional  particles  are  put  into  a  table  along  with 
the  number  N-1.  This  entry  in  the  table  is  defined  as  (N-1)  latent  particles. 

The  tracking  is  continued  with  the  N^^  particle.  When  tracking  of  a  particle  is 
terminated,  the  program  searches  the  latent  table,  picks  the  last  latent  particle 
entered  and  starts  tracking  it.  The  particle  history  is  terminated  when  tracking 
of  the  original  particle  and  all  the  latent  particles  is  completed. 

The  region  in  which  the  particle  makes  a  collision  is  determined  as  follows. 
Starting  from  region  Rj,  the  above  procedure  is  continued  until  it  is  found  that 
in  some  region  Rj^ 

A  -  Mt(E,Rj^)  Sk  <  0  (59) 

The  program  then  proceeds  to  bring  the  particle  to  the  collision  point  in  Rj^.  In¬ 
deed,  the  distance  from  the  edge  of  region  Rj^  to  the  collision  point  is  given  by 
hq'vEjRi^)* 

We  must  now  determine  with  which  element  in  the  region  Rj^  the  particle  collides, 
the  type  of  collision,  and  the  energy  and  direction  of  the  emerging  particle,  if 
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any.  This  is  obtained  by  the  interaction  routines  which  are  described  in  Sec¬ 
tion  7.2, 

The  tracking  of  the  particle  coming  out  of  collision  is  resumed  by  choosing  a  new 
A  and  repeating  the  process  described  above. 

7.2  THE  INTERACTION  ROUTINES 

The  interaction  experienced  by  the  particle  is  determined  in  two  steps.  First  the 
interacting  element  is  chosen;  second,  the  specific  event  is  chosen.  If  the  event 
is  scattering,  the  energy  and  direction  of  the  scattered  particle  are  evaluated. 

The  interacting  element  is  chosen  by  comparing  the  product  of  a  random  number 
with  the  total  macroscopic  cross  section  of  the  region  composition,  ^/i»p(E,R), 
with  the  total  macroscopic  cross  section  of  each  element  in  that  composition, 
Atrp(E,ej),  where  ej  denotes  the  element  in  the  composition  of  region  R, 
j  =  1,  2,  ...,  jC,  ...,  L.  K 

jE-I  £ 

j=l  j=l 

then  the  interacting  element  is 

The  type  of  interaction  is  similarly  determined  by  comparing  a  random  number 
to  the  probabilities  that  a  given  event  takes  place.  These  probabilities  have  been 
tabulated  for  each  element  by  the  GENPRO  routine. 

The  ADONIS  subroutine  which  is  used  to  select  the  interacting  element  and  the 
type  of  interaction  is  called  Data  Retrieval,  entry  3  (DR, 3). 

The  incident  energy,  E,  is  referenced  with  respect  to  the  output  mesh  energies 
Ep  in  an  energy  interval  (Ep,Ej^i)  by  (DR,1).  For  each  element,  at  each  Ep, 
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GENPRO  supplies  tables  of  total  microscopic  cross  sections,  scattering  probabili¬ 
ties  ,  absorption  probabilities,  and  a  scattering  index.  This  index  gives  the  type 
of  scattering,  if  any,  e.g.,  Uastic,  discrete  or  continuum  inelastic.  If  elastic 
scattering  is  indexed,  it  gives  the  type  of  elastic  scattering  that  the  particle  suf¬ 
fers.  For  anisotropic  elastic  scattering,  discrete  inelastic,  or  continuum  inelastic 
scattering,  the  scattering  index  gives  the  proper  location  of  the  tables  of  probabil¬ 
ities  supplied  by  GENPRO:  CHI,  PLEV,  and  ENN  tables  (see  Section  3). 

/ 

(DR,3)  compares  a  random  number,  in  (0,1)  to  the  probability  of  scattering  at 
energy  E  [obtained  by  interpolation  using  the  f  computed  by  (DR,1)]. 

If  the  event  is  absorption,  the  number  of  absorptions  in  the  element  and  the  region 
are  updated,  the  particle  tracking  is  terminated,  and  tracking  of  a  latent  particle, 
if  any,  is  started. 

If  the  event  is  scattering,  the  energy  after  scattering,  EPRIM,  (E'),  and  the  cosine 
of  the  angle  of  scattering  in  the  laboratory  system,  CSTHT,  (cos  0),  are  computed 
by  calling  the  appropriate  routines:  ELAS  and  INELAS  for  neutron  elastic  and 
inelastic  scattering,  respectively,  and  NIKI  for  gamma  elastic  scattering. 

ELAS  Routine 

The  control  digit  NCBD  is  set  to  1,  2,  3,  4  for  the  following  cases:  (1)  isotropic 
in  center  of  mass  system,  (2)  hydrogen  scatterer,  (3)  anisotropic  scattering,  and 
(4)  isotropic  in  laboratory  system  -  no  energy  degradation. 

1.  Isotropic  Scattering  in  Centeii  of  Mass  System 

It  can  be  shown  that  if  ^  is  random  in  (0,1) 

X  =  ^c.,  =  ^  (61) 

where  6^  is  the  scattering  angle  in  the  center  of  mass  system.  E'  and 
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cos  9  are  then  obtained  by  the  standard  relations 


E'  =  E 


(62) 


cos  6  = 


A+1-2AX 


(A+l) 


4AX 

(Sn? 


T/T 


(63) 


where  A  =  the  element  atomic  weight 
E  =  the  incident  energy. 

2.  Hydrogen  Scatterer 

Upon  collision  with  hydrogen,  particles  are  scattered  in  the  foi’ward 
direction  only.  It  can  be  shown  that  the  probability  distribution  function 
cos  6  is  given  by  2  cos  6. 

To  pick  a  cos  6  from  the  distribution  2  cos  6,  two  random  numbers^^  are 
selected  and  the  larger  of  the  two  is  retained  as  cos  6.  Then 


E'  =  E  (cos  ef 


(64) 


3.  Anisotropic  Scattering 


The  cosines  of  angles  between  which  elastic  scattering  is  equally  proba¬ 
ble  for  the  specified  energy  mesh  points,  Ep,  are  supplied  in  the  CHI 
tables  referenced  by  the  scattering  index. 

If  the  incident  neutron  energy,  E,  lies  in  the  interval  (Ep,Ep.j),  the  CHI 
table  relative  to  Ej  is  selected.  This  tabulation  consists  of 


Xn 


1  -  cos  ©n 
2 


for  n  =  0,  1,  2,  ...,  N 


where  Xo  =  0.  Xn  =  ^  <  Xn  '  1- 
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In  ADONIS,  N  is  restricted  to  10.  The  proper  x  is  selected  by  choosing 
a  random  number,  and  comparing  the  product  10^  to  the  numbers 
1,  2,  N  =10  successively. 

Assume  that  n  <  10|  <  n+1 
This  means  tha’  Xn  <  X  <  Xn+i 

The  integral  part  of  the  product  10^  is  rejected  and  the  remainder  is 
taken  as  tJie  interpolation  factor,  f.  [This  assumes  a  uniform  probability 
distribution  of  x  in  each  interval  ( Xjj»  Xn+i^*  ] 

Then 

X  =  Xn  +  «X„+,  -  Xn)  (65) 

This  function  is  performed  by  the  Data  Retrieval,  entry  6  routine  (DR,  6). 
Therefore, 

cos  0  =  1  -  2  X  (66) 

and  E'  is  obtained  using  Eq.  82. 

4.  Isotropic  in  Laboratory  System  -  No  Degradation 

If  the  scatterer  atomic  weight  is  large,  Eqs.  62  and  63  show  that 

E'  =  E  (67) 

cos  6  » 1  -  2x  =  cos  (68) 

Thus,  the  scattering  is  approximately  isotropic  in  the  laboratory  system 
and  the  particle  energy  remains  unchanged  by  the  collision.  Therefore, 
a  random  number,  is  chosen  and  cos  0  =  1-2^. 
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INELAS  Routine 


Data  Retrieval,  entries  7  and  8  subroutines,  (DR, 7)  and  (DR, 8)  are  entered  for  dis¬ 
crete  and  continuum  scattering,  respectively  (overlapping  is  included  in  the  con¬ 
tinuum  spectrum),  to  select  E'. 

If  E  lies  in  (Ep,Ej^j),  the  PLEV  and  ELEV  tables*  for  discrete  scattering  or  the 
ENN  table*  for  continuum  scattering  relative  to  Ep  are  selected. 

To  select  from  the  ENN  table,  (DR, 8)  uses  the  same  technique  as  (DR,6). 

To  select  from  the  PLEV  table,  (DR, 7)  uses  the  same  technique  as  (DR, 3)  in  the 
selection  of  the  interacting  element. 

Using  the  inequality  (Expression  60)  in  which  the  element  e^  is  replaced  by  the 
energy  levels,  one  obtains  the  appropriate  energy  level  to  which  the  particle 
scatters  and  E'  =  E  -  E*^. 

The  NIKI  Routine 

It  is  assumed  that  y-scattering  results  from  the  Compton  effect.  A  y-ray  of 
energy  a  (in  mc^)  traversing  matter  has  a  probability,  p(x/or)  dxdf,  of  having  a 
collision  in  the  distance  df  and  emerging  from  this  collision  with  an  energy 
a'  =  xof.  For  x  satisfying  1  <  x  20, 

P(x/a)  dx  =  :4  (cos'^  0  -  1  +  X  +  ^  dx.  (69) 

where  A  =  mal/ a 

n  =  the  number  of  electrons  per  cm^ 

00=  the  classical  radius  of  the  electron 

9  =  the  scattering  angle  which  is  related  to  o  and  x  by  the  Compton  relation 
■•See  Section  3. 
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(70) 


1  y 

cos  0  =  -  -  -  +  1. 
a  a 

p(x/a)  is  called  the  Klein-Nishina  probability  distribution  function.  To  pick  a  x 
from  p(x/«)>  a  rejection  technique  devised  by  H.  Kahn  is  used  in  ADONIS.*^ 
pix/ot)  is  broken  up  into  the  sum  of  two  probability  distribution  functions.  A  ran¬ 
dom  number,  is  compared  to  (2a+l)/(2Qr+9). 

If  <  (2a'+l)/(2Q!+9),  y  is  set  equal  to  1  +  Zaii  (where  ^2  is  a  second  random 
number),  and  another  random  number,  is  compared  to  4  [(1/y)  “(1/y^)]  =  na. 

If  it  is  less  than  m,  x  =  y*  K  it  is  greater  than  m,  y  is  rejected  and  the  procedure 
is  repeated. 

If  >  (2a+l)/(2Q;+9),  y  is  set  equal  to  (2Q!+l)/(l+2a^2)»  ^3  t®  compared  to 

l(cos«ea), 

where 

1  V 

cos  6  =  —  -  +  1. 
a  a 

If  it  is  less  than  (1/2)  [cos*  o  (1+y)],  x  =  y*  If  it  is  greater  than  (l/2)[cos*  9+ 
(Vy)]»  y  is  r^ected  and  the  procedure  is  repeated. 

Once  X  has  been  selected,  NIKI  proceeds  to  compute  and  cos  6. 

The  interaction  routines  compute  and  transmit  to  the  main  program  the  energy 
after  scattering,  E',  and  the  cosine  of  the  scattering  angle,  cos  6.  What  is  re¬ 
quired  is  to  evaluate  the  direction  cosines  of  the  new  particle  direction  in  the 
Cartesian  coordinate  system,  given  cos  6  and  the  initial  direction  cosines.  The 
DIREC  routine  computes  these  angle  cosines  by  means  of  two  subroutines; 


1.  Subroutine  PHI  uses  a  rejection  technique  to  choose  the  cosine  and  sine  of 
an  azimuthal  angle  equidistnbuted  in  the  interval  0  <  0  <  27r:  cos  0  and 
sin  0. 

2.  Subroutine  DCB  computes  the  direction  cosines 

fix  =  ^  sin  0  )  +  fixM 

fiy  =  sin  0  -  fix  cos  0)  +  fiyM 

fi^=  /(iV)  (l-fi|)  sin0  +  fizM  (71) 

where  m  =  cos  ^ 

is  the  given  initial  direction. 

7.3  VARIANCE  REDUCTION  TECHNIQUES 

As  pointed  out  in  the  introduction  to  this  section,  the  transport  game  is  altered 
for  better  efficiency.  In  shielding  problems,  one  is  interested  in  particles  from 
a  source  which  are  transmitted  through  the  shield.  In  problems  of  interest,  the 
transmission  is  generally  small,  and,  if  results  with  reasonable  accuracy  are  re¬ 
quired,  an  extremely  large  number  of  particle  histories  must  be  considered.  This 
involves  large  amounts  of  computing  time;  thus  it  is  costly. 

To  obtain  a  practical  calculation,  one  must  do  other  than  play  the  direct  analog 
game  in  which  a  great  deal  of  time  is  spent  on  particles  which  are  not  trans¬ 
mitted.  Many  methods  have  been  tried  in  which  the  situation  is  improved  by  com¬ 
puting  the  same  average  scores  but  with  a  reduced  variance.  Most  of  *,hese  in¬ 
crease  the  work  per  game,  but  hopefully  result  in  a  net  saving  of  work.  The  meth¬ 
od  of  reducing  variance  in  the  main  program  of  ADONIS  is  a  particular  case  of 
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importance  sampling.  However,  for  the  special  problem  of  evaluating  neutron 
fluxes  in  a  small  volume  detector,  a  different  approach  is  used.  Here,  instead  of 
decreasing  the  variance  by  altering  the  game,  the  variance  is  decreased  by  chang 
ing  the  method  of  scoring. This  method  is  called  statistical  estimations,  or 
scoring  by  analytical  estimations. 


7.3.1  ADONIS  Importance  Sampling 

In  altering  the  direct  analog  game  of  chance  to  reduce  variance,  one  tries  to  esti¬ 
mate  the  desired  result  by  emphasizing  the  choice  of  phase  space  points  which  lead 
to  higher  expected  contribution  to  the  score.  I(x,E,fi)  is  a  measure  of  the  “im¬ 
portance”  of  the  point  (2c,E,fl). 


It  can  be  shown  that^  the  use  of  I(3^E,fD  leads  to  a  game  in  which  particles  are 
picked  from  an  altered  source,  S(2t,E,fy,  and  collision  locations  and  energy  and 
direction  of  scattered  particles  are  chosen  from  an  altered  kernel,  K(x— x';E— E'; 
such  that 


5/v  F 

/  I(2t,E,^  S(x,E,iy  dx  dE  dn 


and 


(72) 


k(x-x';E.-E';fi-.i2') 


(73) 


where  S  and  K  are  the  normalized  source  and  the  collision  kernel,  respectively, 
of  the  integral  Boltzmann  equation. 


The  questions  are  how  to  choose  the  most  efficient  importance  function,  I,  and, 
having  chosen  I,  how  do  we  play  the  game  of  chance? 

A  proper  choice  of  I  forces  particles  to  points  in  phase  space  such  that  their  con¬ 
tribution  to  the  desired  answer  will  be  large.  As  a  matter  of  fact,  there  exists 
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an  ideal  importance  function  which  leads  to  a  zero  variance  game.  It  can  be  shown 
that  the  most  efficient  important  functions  are  solutions  to  the  adjoint  equations 
to  the  Boltzmann  equation. Obtaining  these  solutions  presents  an  equally  dif¬ 
ficult  problem.  In  practice,  the  choice  of  I  is  limited  because  it  depends  on  the 
complexity  of  the  kernel,  K.  Tue  Boltzmann  equation  kernel,  K,  is  itself  the 
product  of  the  transport  kernel,  T(x-»x';  E,  fi) ,  and  the  collision  kernel, 

C(x';  E-E';  Hence,  we  could  attach  an  importance  function  to  each  of 

these  kernels  separately.  The  “importance”  of  a  particle  will  then  be  governed 
by  the  importance  of  the  region  in  the  phase  space  in  which  it  finds  itself  at  the 
moment.  In  ADONIS,  the  importance  function  is  taken  as  a  piecewise  constant 
function  of  energy  and  space  and  independent  of  direction  angle,  i.e.,  I(E,x)  = 
constant  in  each  region  Ri  and  each  energy  interval,  Ej  <  E  <  Ej^j. 

Consider  now  a  particle  at  energy  E  leaving  collision  at  x  in  region  R^.  When 
the  particle  goes  from  x  to  x',  the  collision  density  is  changed  by  the  ratio  of  the 
importance  function  at  the  new  point  to  the  importance  function  at  the  old  point, 
[l(E,x')]/[l(E,x)].  To  insure  the  correct  collision  density  at  the  new  point,  we 
may  follow  rhe  particle  from  x,  to  x'  using  T(3c— x';E).  Having  reached  the  new 
collision  point,  create  the  number  [l(E,x')]/[l(E,x)]  to  go  into  collision  at  jc'.  T  is 
method  assures  us  ol  always  having  the  correct  collision  density  since 

T(x-x';E,fl)=  T(x^x';E,n)  (74) 

In  ADONIS,  a  slightly  different  method  is  used  by  considering  the  region  bounda¬ 
ries  at  which  the  importance  function  changes.  Suppose  that  the  particle  starting 
from  X  in  Ri  crosses  the  boundary  between  regions  Rj  and  Rj+j  and  that 

I(E,Ri) 

Then  the  collision  density  is  supposed  to  double.  One  way  of  insuring  this  is  by 
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creating  another  particle,  precisely  like  the  first,  and  then  following  both.  Of 
course,  when  the  particle  crosses  the  boundary  from  Ri+j  to  Rj[,  we  have  to  halve  it. 
One  way  of  halving  a  particle  is  by  “killing”  it  half  the  time.  The  first  process  is 
called  “particle  splitting”  and  the  second  process  is  referred  to  as  “Russian 
roulette.  ’  ’ 

Similarly,  when  the  particle  enters  collision  at  x',  say  in  region  Rj^,  the  “splitting” 
or  the  “Russian  roulette”  technique  is  used  according  to  the  magnitude  of  the 
ratio  I(E',Rjj.)/l(E,Rjj).  We  could  make  the  importance  function  resemble  even 
more  closely  the  true  function  by  including  an  angular  dependence.  However,  it 
is  no  small  task  to  assign,  even  foi  simple  geometries,  an  importance  function 
as  a  function  of  ^  E,  and  in  each  region. 

Note  that  a  weighting  function  is  used  in  ADONIS.  The  weight  is  merely  the  in¬ 
verse  of  the  importance,  i.e., 

™=T(I;r) 

When  a  particle  goes  from  a  point  [E,x(Rj)]  to  a  point  [E^x(Rj)] ,  we  create  or  kill 


particles.  The  technique  used  in  ADONIS  to  create  or  kill  particles  has  been 
described  in  Section  7.1. 

We  now  investigate  techniques  whereby  one  can  get  relatively  good  importance 
functions. 

For  the  situation  in  which  the  importance  function  is  dependent  on  position  only, 
I(R),  it  has  been  shown  that  it  shou  d  have  the  property  of  keeping  the  particle 
flux  about  constant  from  one  region  to  the  next.  What  is  done  in  practice  is  to 
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run  a  problem  with  relatively  few  histories  and  see  whether  the  condition  is  satis¬ 
fied.  If  it  is  not,  then  we  alter  the  boundaries  according  to  the  following  criterion. 
Whenever  the  number  of  particles  falls  off  too  rapidly,  one  introduces  additional 
importance  boundaries,  ox  increases  the  value  of  the  importance  functions  and 
keeps  the  boundaries  intact.  The  reverse  is  done  whenever  the  number  of  particles 
increases  too  rapidly.  Some  trial  and  error  along  with  educated  guessing  can 
often  lead  to  a  good  importance  function. 

When  the  importance  function  depends  upon  energy  as  well,  the  same  method  is 
used  by  recording  the  number  of  particles  crv^ated  in  ep.cn  energy  group  for  each 
spatial  region. 

7.3.2  Scoring  by  Analytical  Estimations  -  Flux  at  a  Point^* 

In  some  neutron  calculations,  it  is  desirable  to  obtain  answers  for  flux,  dose, 
or  the  like  at  a  point.  Since  no  history  can  be  expected  to  carry  a  neutron  through 
a  point,  it  is  necessary  to  use  anal5rtic  expected  value  estimates.  The  method 
uses  the  relation  between  E, dx'  dE  df^  the  density  of  particles  er.'tering 
collision  in  volume  (k'  with  energy  in  dE  and  direction  in  dfi^  and  0(i),  the  flux 
at  x. 

'  /  ,  \  [  -u(E')  ix-x'l]  ' 

^.(x)  =  /  dE'  ;  (tc'  /  dE  /  dn  (Ki',E,!3  X  j2g  •  njEj  J  ■ 

(71) 

where  g(w;E)  dw  =  the  probability  of  scattering  of  a  neutron  of  energy  E  through 

an  angle  whose  cosine  is  in  the  range  (w,  w+dw) 

E'  =  the  energy  after  scattering 

M(E')  =  the  total  attenuation  coefficient  of  the  scattered  neutron. 
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The  problem  is  to  evaluate  the  estimator  in  brackets  in  Eq.  77  every  time  a  par¬ 
ticle  goes  into  collision  and  average  it  over  all  histories.  A  straightforward 
procedure  yields  an  infinite  variance  if  Ix-x'l  can  become  very  small.  Indeed, 
the  estimator  has  a  ix-x'  1“^  behavior.  The  method  devised  by  Kalos  and  em¬ 
ployed  in  ADONIS  to  avoid  infinite  variances  is  to  remove  the  singularity  from 
the  scoring  function  and  to  incorporate  it  in  the  sampling  probability  distribution 
function  (pdf). 

What  is  required  is  that  for  every  collision,  we  estimate  the  “once-more  collided” 
contribution  to  the  flux  at  a  point  and  choose  the  position  of  the  intermediate  col¬ 
lision  to  make  the  variance  finite.  Hence,  the  final  estimate  is  the  collided  flux 
only.  A  separate  analytic  calculation  must  be  made  for  the  uncollided  flux. 

H  S,  D,  and  A  denote  the  collision  point,  the  detector,  and  the  intermediate  point, 
respectively,  and  R  =  ISDl,  Ti  =  ISA  I,  and  r2  =  lADN  the  following  pdf  is  select¬ 
ed: 


UjdV  = 


8n 


-M  rj  -ii'r2 


dV 


(78) 


where 


M'=M(E)+i  (79) 

Noting  that  U2  is  an  average  of  two  pdf’s,  each  of  the  foi  m 

/  -M'r 
M  e 

4n  r^ 

where  the  origin  of  r  is  taken  either  at  the  collision  point  or  the  detector  point, 
it  is  easy  to  sample  from  U2. 
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1.  Choose  with  equal  probability  either  S  or  D 

2.  Choose  the  direction  toward  A  isotropically 

•  i/r 

3.  Choose  a  distance  r  from  the  pdf  p'  e 

This  distance  is  either  r^  or  r2  and  the  remaining  distance  must  be  computed. 
(Note  that  the  direction  toward  the  intermediate  point  A  is  chosen  from  an  iso¬ 
tropic  distribution.  Thus,  if  the  source  particle  angular  distribution  is  not  iso¬ 
tropic,  the  first  direction  picked  from  S  must  satisfy  this  distribution.  At  the 
present  time,  ADONIS  assumes  isotropic  angular  distribution  of  the  source  par¬ 
ticles  when  the  flux  at  a  point  routine  is  used.) 

The  contribution  to  the  “once  more  collided  flux  at  a  given  detector  D’’  from  a 
neutron  of  energy  E  and  direction  entering  collision  at  S  is  then  given  by 

-S  p(Ei,Rj)Si 

0(E{,D)  =  {[P(E,S)]  g(g-I?i;E)  e 


-2:/i(E2,Rj)Sj 

[P(Ei,A)]MTfEi,R(A)]  e  ^  } 

/U2(4nrj)(4nrl)  (80) 

where  i2i,  Ej  =  the  direction  and  energy  after  scattering  at  S 
E2  =  the  direction  and  energy  after  scattering  at  A. 

r  p  r, 

where  S^,  i  =  1,  2,  ...,  denote  the  track  lengths  in  the  regions  that  the  particle 
traverses  from  S  to  A. 


S  Si  =  r2 

i  ’ 


IIS 


where  Sj,  j  =  1,  2,  denote  the  track  lengths  in  regions  between  A  and  D.  U2 
is  computed  using  Eq.  78. 

[P(E,S)]  means  that  a  game  of  chance  is  played  to  find  the  interacting  element  and 
the  type  of  interaction. 

K  the  interaction  at  S  or  A  is  an  absorption,  the  contribution  to  the  once  more 
collided  flux  at  D  is  zero.  If  it  is  scattering,  Eq.  80  is  computed. 

The  contribution  to  the  uncollided  flux  is  finally  added  to  obtain  the  total  flux  at 
detector  D  in  each  output  energy  interval. 

“]^  M(E2,Rk)  Sjj 


where 


S  Sk  =  R 

k 

where  k  =  1,  2,  ...,  denote  the  track  lengths  in  the  regions  between  S  and  D, 
and  E2  means  tiiat  E2  is  in  the  output  energy  interval. 

7.4  TYPES  OF  ANSWERS 
The  ADONIS  output  consists  of: 

1.  The  average  flux  per  region  and  per  unit  energy  in  each  output  energy 
interval,  and  its  deviation 

2.  The  volume  V,  the  weight  W(E,R),  the  number  of  births  (latent  particles), 
deaths  (by  importance  sampling),  energy  degradations  (particles  whose 
energy  falls  below  a  specified  cutoff  energy,  if  any)  and  absorptions  in 
each  region  except  the  transmitting  and  reflecting  regions 

3.  The  number  of  transmitted  and  reflected  particles 
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4.  As  an  option,  a  total  fine  response  and  its  deviations 

5.  As  an  option,  co^,rse  responses  for  any  region 

6.  An  interaction  tape  for  neutron  problems  only 

7.  A  transmitted  particle  tape. 

Items  2,  3,  6,  7  have  been  discussed  in  the  above  subsections.  This  description 
will  be  restricted  to  items  1,  4,  and  5. 

7.4.1  Average  Flux 

There  are  two  ways  to  measure  the  average  flux,  in  a  region.  The  usual  way 
is  to  count  m(R),  the  number  of  interactions  in  the  region,  and  if  N  is  the  number 
of  source  particles. 


0  = 


m(R) 

mt(R)  n 


(82) 


This  method  is  inefficient  when  R  is  a  thin  region.  It  is  then  necessary  to  use 
the  other  method  which  is  called  the  track  length  method: 

-  _  Total  distance  all  particles  traveled  in  R 

^  (Volume  of  R)  X  N  ^ 

It  can  be  shown  that  both  methods  are  equivalent.  ADONIS  computes  average  flux 
per  region  and  per  unit  energy  in  each  output  energy  interval  using  the  track 
length  method. 

The  source  particles  are  considered  in  G  equal  groups  of  n  particles  each,  such 
that  N  =  Gn,  each  group  being  denoted  by  g  =  1,  2,  ...,  G.  For  most  cases, 

20  particles  per  group  are  taken.  In  any  case,  n  must  be  specified  in  the  input. 
This  grouping  is  made  for  statistical  purposes.  As  it  will  be  shown,  although 
the  contribution  to  the  average  flux  of  all  track  lengths  is  cumulative,  the  con¬ 
tribution  to  the  variance  depends  strongly  on  the  grouping. 
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Let  Sg(R,E^)  be  the  sum  of  all  track  lengths  of  particles  from  the  group  in 
region  R  with  energy  E  in  the  output  energy  interval. 

The  average  track  length  per  soirrce  particle  per  unit  energy  in  the  interval 
is  given  by 


S(R,Ej)  = 


G 

2  Sg(R,Ej) 
g=^ 

xGn 


(84) 


where  Ej,  Ej^j  defir.  the  energy  interval. 

The  average  flux  per  source  particle  per  unit  energy  in  the  interval  in  region  R 
is  given  by 


Fh  -  S(R>EP)  W(Ej,R) 

>  -  v(R) 


(85) 


where  V(R)  and  W(EJ,R)  are  the  volume  and  the  weight  in  tlie  energy  group  of 
region  R,  respectively.  If  the  energies  are  given  in  ev,  track  lengths  in  cm,  and 
V(R)  in  cm®,  ^(E^,R)  is  obtained  in  number  of  particles  per  ev-cm®-source  par¬ 
ticle. 


If  V(R)  is  infinite,  the  ADONIS  output  is  the  product  S(R,Ej)  W(Ej,R). 


The  deviation  of  the  average  flux  is  obtained  from 


(7(R,Ej)  = 


0(Ej,G)»G 

G 

S 

g=l 


G 

G 

2' 

2  Sg(R,Ej) 

J 

g=l 

G 

G 

j 

1/2 


(86) 
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7.4.2  Fine  Response 


If  a  functional  of  the  type 

f(H)  =  /r  /e  f(}^E)  dX  dE  (87) 

is  desired  (fine  dose,  heating  ...),  the  ADONIS  input  contains  a  table  of  f(R,Ej  ) 
at  the  energies,  Ej^,  of  the  cross  section  used  in  ADONIS  (81  energy  mesh  points). 
Because  of  core  storage  restrictions,  a  maximum  of  five  possible  distinct  re¬ 
gions  can  be  considered  for  fine  response  evaluation  (5  x  81  words  for  f^R,E)  are 
available]. 

Each  particle  track  length,  T(R,E),  is  multiplied  by  an  interpolated  value  of  the 
tabulated  fine  response  function,  f(R,E),  and  the  products  are  accumulated.  Then 
the  total  fine  response,  f(R)  is  given  by 


f(R)  = 


W(R) 

G-n-V(R) 


2  T(R,E)  f(R,E) 
g=l 


The  deviation  of  the  fine  response  is  given  Dy 

/  f  G  f  G 

,  W(R)  /  J_  I  _ gfl _ 

^  n-V(R)\G-l  I  G  G 


21  \  1/2 


Note  that  at  the  present  time,  a  fine  response  calculation  can  be  performed  only 
if  the  region  weight  is  not  energy  dependent. 


7.4.3  Coarse  Response 

In  any  region,  f  can  be  obtained  by  giving  as  input  tables  of  f(E^,R),  where 
f(E^,R)  is  piecewise  constant  in  each  output  energy  interval.  The  coarse  re- 
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sponse  is  then  given  by 


f  = 


0(R,Ej)  f(Ei,R) 


(90) 


7.5  INTERMEDIATE  RESULTS 

For  debugging  purposes  and  other  reasons,  a  trace  using  sense  switch  3  has  been 
built  into  the  ADONIS  code.  This  trace  permits  one  to  follow  the  flight  of  a  given 
number  of  particles.  When  the  number  of  specified  source  particles  is  exhausted, 
the  program  calls  the  edit  routine  STAT. 


ADONIS  problems  may  run  for  long  periods  of  time.  In  order  to  save  problems  so 
that  they  may  be  restarted,  the  parameters  used  by  ADONIS  have  all  been  stored 
in  COMMON,  so  that  if  COMMON  were  written  on  a  tape,  read  back,  and  the  pro¬ 
gram  entered  at  the  proper  place,  one  could  continue  a  problem.  This  coding 
has  been  placed  in  the  ADONIS  code  and  is  activated  by  sense  switch  1  whose 
function  is  to  dump  COMMON  and  continue  the  problem  (the  computer  will  print 
a  message  telling  the  operator  to  turn  off  sense  switch  1  and  will  hang  waiting 
till  he  does).  Sense  switch  1  and  sense  switch  6  set  together  will  dump  a  problem, 
force  an  edit,  and  terminate.  To  pick  up  a  problem  that  has  been  dumped,  one 
places  a  1  in  the  restart  option  of  the  first  parameter  card  and  resubmits  the  job 
as  a  normal  ADONIS  run.  The  code  reads  the  dump  tape  and  picks  up  from  where 
it  left  off.  The  first  parameter  card  must  be  exactly  the  same  as  the  initial  param¬ 
eter  card  with  the  following  two  exceptions:  (1)  the  last  history  number  may  be 
increased  and  (2)  the  restart  option  must  be  a  1.  The  title  and  the  first  card  are 
all  the  input  necessary  for  restarting  a  problem.  If  coarse  response  function  data 
are  present,  these  must  follow  the  first  parameter  card. 
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7.6  ADONIS  INPUT 


The  input  required  for  ADONIS  consists  of  source  input,  cross-section  input, 
geometry  input,  and  parameter  input  describing  the  particular  problem  of  interest. 

Source  Input 

The  source  input  consists  of  a  source  tape  containing  source  particles  generated 
by  programs  like  VANGEN  or  GASP.  These  source  particles  are  arranged  by 
particles  blocked  in  groups  of  100;  each  particle  is  represented  by  its  history 
number,  its  initial  position,  direction,  and  energy,  and  the  region  of  birth. 

Cross-Section  Input 

The  cross-section  input  consists  of  the  desired  cross  sections  arranged  in  tables 
so  that  they  are  accessible  to  the  ADONIS  interaction  routines.  These  tables  are 
placed  in  the  COMMON  portion  of  memory  by  the  program  DATORG.  This 
COMMON  is  consistent  with  that  used  by  the  ADONIS  interaction  routines.  A 
variable  storage  design  has  been  employed  so  that  only  that  storage  necessary 
for  a  given  problem  is  reserved,  leaving  the  remaining  storage  available  for  use 
where  necessary  so  that  the  originator  may  be  able  to  expand  the  input  where  de¬ 
sired. 

Geometry  Input 

Tables  1,  2,  3  described  in  EZGEOM  are  transmitted  to  ADONIS  via  COMMON. 
Problem  Dependent  Parameter  Input 

The  program  MPINP  reads  in  cards  describing  the  problem. 

The  first  parameter  card  contains  the  limits  of  the  problem  as  shown  below. 

Entry  Type  of  Limit 

1  The  number  of  the  first  history  to  be  treated  (usually  1) 

2  The  number  of  the  last  history  to  be  treated 
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Entry 


Type  of  Limit 


3  The  number  of  flux  bins  for  flux  tracking 

4  The  number  of  particles  per  group  for  statistical  scoring  of 

flux 

5  The  total  number  of  ordinary  and  nonordinary  regions  in  the 

problem 

6  The  transmitting  region  number 

7  The  total  number  of  distinct  weights,  Q 

8  The  number  of  distinct  region  dependent  weights,  q 

9  The  number  of  distinct  compositions  (given  as  input  to  the 

associated  DATORG  problem) 

10  The  number  of  flux-at-a-point  detectors  (^15),  if  any 

11  The  number  of  fine  response  functions  ( <6) 

12  Restart  option: 

0  =  original  problem 
1  =  restart 

13  The  type  of  problem  to  be  run: 

0  =  neutron 
1  =  gamma 

Cards  following  the  first  parameter  card  contain  the  items  listed  below. 

1.  The  cutoff  energy  and  the  thermal  energy,  if  any.  If  the  thermal  energy 
entry  is  zero,  particles  slowing  down  below  the  energy  cutoff  are  termi¬ 
nated  and  counted  as  degraded  particles.  If  the  thermal  energy  entry  is 
nonzero,  particles  with  energy  below  the  energy  cutoff  are  not  terminat¬ 
ed.  They  are  assumed  to  be  at  thermal  energy  and  no  longer  lose  ener¬ 
gy  upon  scattering.  The  cross  sections  used  are  the  thermal  cross  sec¬ 
tions.  Absorption,  reflection,  transmission,  or  death  by  importance  are 
then  the  only  ways  a  history  may  be  terminated. 

2.  The  energy  mesh  points  which  define  the  output  energy  intervals  follow. 
The  number  of  entries  corresponds  to  the  number  of  specified  energy 
bins  plus  one.  The  entries  are  in  ev  and  in  descending  order. 

3.  The  next  input  entries  are  the  atomic  weights  per  composition.  This  in¬ 
put  consists  of  the  atomic  weights  for  each  element  of  each  specified  com¬ 
position,  one  composition  per  card.  These  cards  are  arranged  in  the 
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order  corresponding  to  the  DATORG  input.  This  input  is  used  as  a 
check  on  ihe  cross-section  input  read  from  the  DATORG  output.  When 
a  particle  collides  and  the  interacting  element  is  determined,  the  mass 
of  the  element  is  taken  from  this  list. 

4.  Next  the  weight  numbers  which  are  assigned  to  each  region  appear. 

Region -dependent  weights  and/or  sets  of  region-  and  energy -dependent 
weights  must  be  assigned  to  each  region  in  the  problem,  as  shown  in 
Fig.  7.  (If  no  importance  sampling  is  desired,  1.0  is  assigned  as  the 
weight  for  all  regions.)  For  each  distinct  region-dependent  weight  and/ 
or  set  of  energy-  and  region-dependent  weights,  a  unique  number  is  as¬ 
signed  as  follows:  first,  the  numbers  1  to  q,  where  q  is  the  number  of 
distinct  region -dependent  weights,  are  assigned  to  each  of  these  weights; 
then  the  numbers  q+1  to  Q,  where  Q  is  the  total  number  of  distinct  weights, 
are  assigned  to  each  of  the  sets  of  energy-  and  region -dependent  weights. 
This  procedure  is  shown  on  the  table  in  Fig.  7. 

The  first  set  of  cards  lists,  in  the  order  of  the  region  number  (14  to  a 
card),  the  weight  number  assigned  to  each  region,  as  shown  in  Fig.  7. 

In  the  next  set  of  cards  the  q  region  weights  are  entered  in  the  order 
that  the  weight  numbers  w’ere  assigned.  Then  the  Q-q  sets  of  energy - 
dependent  weights  are  entered,  each  set  starting  a  new  card  and  in  the 
order  of  their  weight  numbers  (Fig.  7). 

5.  The  composition  number  per  region  input  follows.  The  input  consists 

of  the  composition  number  from  one  to  the  number  of  different  composi¬ 
tions  in  the  problem  being  assigned  to  each  region.  (This  is  the  same 
procedure  as  that  described  for  the  weight  numbers.) 

6.  The  next  type  of  input  is  optional.  The  number  of  fine  response  functions 
has  been  specified  in  the  first  parameter  card.  Each  response  function 
contains  a  title  card,  a  card  denoting  the  region  to  which  this  response 
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function  applies,  and  the  response  function  enumerated  at  each  of  the 
mesh  point  energies,  Ep,  that  were  used  to  create  the  cross -section 
tape. 

7.  The  following  input  for  flux  at  a  point  is  optional.  This  input  consists  of 
a  card  containing  the  number  of  output  energy  bins  plus  one.  The  limits 
of  the  energy  intervals  are  listed  in  ev  and  in  descending  order.  A  title 
card  is  next.  Following  the  title  card  are  cards  containing  the  x-y-^  co¬ 
ordinates  of  the  detector  points  and  the  regions  in  which  the  detectors 
are  located.  The  number  of  these  cards  corresponds  to  the  number  of 
detectors  specified. 

8.  Finally,  the  coarse  response  function  input  appears.  This  input  consists 
of: 

a.  A  card  denoting  the  number  of  coarse  response  functions  desired 
(0  means  no  coarse  response  is  desired) 

b.  Each  response  function  is  described  in  the  following  cards. 

(1)  A  title  card 

(2)  A  region  information  card  containing  the  number  of  regions  for 
which  this  coarse  response  function  applies 

(3)  The  average  response  function  for  each  output  energy  interval 

(4)  A  table  of  the  regions  to  which  this  response  function  applies. 

7.7  ADONIS  INPUT  FORMS 

ADONIS  input  forms  r.re  shown  on  the  follov*'ing  pages. 
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ADONIS  GROUP-2  INPUT  (CONTINUED) 


ADONIS  GROUP-2  INPUT  (CONTINUED) 


*N<  .  soquciil 


UEGION  DEPENDENT  WEIGHT  CARDS 


(See  F»g.  7,  second  card) 


ADONIS  GROUP-2  INPUT  (CONTINUED) 


(See  Fig.  7,  third  and  fourth  cardsi 


ADONIS  GROUP-2  INPUT  (CONTINUED) 


ADONIS  GROUP-2  INPUT  (CONTINUED) 


TITLE  CARD  FOR  DETECTORS 


^Number  sequentially 
Note: 

These  cards  must  be  In  direct  correspondence  with  energy  bln  cards  In  this  section. 
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8.  DESCRIPTION  OF  THE  UNC-SAM  PROGRAM  SEQUENCER  (SAM) 


UNC-SAM  contains  a  program  sequence*'  '!)  whose  functions  are  listed  below. 

1.  To  read  the  program  name  ca.  ...  appears  in  the  input.  For  example, 
in  Fig.  8,  it  will  read  first  the  DATORG  card. 

2.  To  check  that  programs  that  supply  input  have  been  run.  h  the  case  of 
DATORG,  it  will  check  that  the  output  of  GENPRO  is  available. 

3.  To  read  in  and  execute  the  appropriate  program. 

4.  Upon  completion  of  the  program  to  check  that  no  error  has  occurred.  In 
the  event  of  an  error,  UNC-SAM  will  bypass  all  subsequent  programs 
which  require  the  output  of  this  program. 

5.  To  read  the  next  prc^ram  name  card  and  to  repeat  the  procedure  from 
item  2. 

When  SAM  reads  an  “END”  as  a  program  name  card,  the  program  sequence  is 
terminated.  SAM  is  the  main  program  on  an  overlay  tape  (logical  tape  3)  contain¬ 
ing  the  DATORG,  EZGEOM,  VANGEN,  GASP,  and  ADONIS  programs. 

Input*  required  for  SAM  are  the  numbers  0,  1,  or  2  corresponding  to  each  of  the 
programs  on  the  overlay  tape. 

0  means  the  corresponding  program  will  not  be  run  in  this  sequence. 

1  means  the  corresponding  program  will  be  run. 

♦The  description  of  the  proper  program  input  appears  in  the  input  section. 
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Function 

Response 

oarse 


Fig.  8  —  Input  sequence  and  tape  assignment  for  running  DATORG, 
EZGEOM,  VANGEN,  and  ADONIS 
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2  means  the  corresponding  program  has  been  run,  and  that  the.  output  is 
available  for  a  program  that  will  be  run  in  this  sequence. 

The  GENDA  and  GENPRO  programs  that  supply  the  Element  Data  Tapes  (i.e.,  the 
cross-section  library)  are  separate.  The  Element  Data  Tape  (output  from  GENPRO) 
is  required  to  run  DATORG.  Consequently,  the  numbers  0  or  2  only  apply  as  an 
entry  in  the  SAM  input  for  the  program  GENPRO. 

A  program  in  the  sequence  may  appear  in  any  order,  provided  all  the  programs 
that  supply  input  io  it  have  previously  been  run. 

The  input  required  for  each  prc^ram  is  preceded  by  a  program  name  card. 

The  storage  has  been  arranged  so  that  all  data  that  are  problem -dependent,  or  are 
necessary  for  restarting  a  problem  from  a  previous  dump,  appear  in  the  COMMON 
portion  of  core.  With  this  arrangement  one  can  dump  only  the  COMMON  portion 
of  the  core  on  tape  and  be  able  to  restart  an  ADONIS  problem  (Fig.  9).  The 
COMMON  storage  tape  (logical  9)  must  be  present  at  all  times  when  UNC-SAM  is 
being  used. 

The  program  DATORG  requires  as  input  a  tape  containing  the  library  of  cross 
sections  (neutron  or  gamma  depending  upon  the  problem  specified)  and  composi¬ 
tion  input  cards. 

To  signify  that  the  cross-section  tape  (logical  tape  13)  is  available,  a  2  is  placed 
in  the  SAM  input  card  corresponding  to  the  GENPRO  prc^ram.  To  signify  that 
DATORG  is  to  be  run,  a  1  is  placed  in  the  DATORG  position  of  the  SAM  card.  The 
cross-section  output  of  a  successfully  completed  DATORG  problem  is  placed  on 
the  COMMON  storage  tape  (logical  tape  9). 

The  COMMON  storage  tape  (logical  tape  9)  is  required  by  the  program  EZGEOM 
so  that  the  output  of  a  successful  run  may  be  added  (if  DATORG  has  been  run  pre- 
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viously),  or  initialized  (if  DATORG  has  not  been  run  previously).  Geometry  card 
input  data  are  also  required  to  run  EZGEOM  (see  Section  6). 


The  program  VANGEN  requires  logical  tape  11  as  an  output  tape  (source  tape), 
and  card  input  data  (defining  the  source  distribution)  (see  Section  5.6.1). 

The  program  GASP  requires  that  logical  tape  1 1  be  assigned  as  output  and  logical 
tape  8  from  a  previously  run  primary  problem  be  assigned  as  input.  Also,  card 
input  data  are  required  (see  Section  5.6.2). 

The  program  ADONIS  requires  a  source  tape  (logical  tape  11)  and  the  COMMON 
storage  t?pe  (logical  tape  9)  as  input,  a  transmitted  particle  tape  (logical  tape  14) 
and  an  interaction  tape  (logical  tape  8  for  neutron  problems  only)  as  output. 

Card  input  data  are  also  required  (see  Section  7). 

An  array  called  TAB  3,  dimensioned  for  15,000  words,  has  been  set  aside  to  store 
the  geometry,  cross  section,  and  statistical  scoring  arrays  for  each  region  (births, 
deaths,  absorptions,  degradations,  sum  of  track-lengths  for  each  energj^  bin  in 
each  region,  the  cumulative  flux,  and  the  cumulative  square  flux,  the  weights  and 
composition  numbers). 

This  array  is  arranged  in  the  following  manner; 


Type  of  Data 


Number  of  Computer  Words  Used 


1.  Geometry  data 

2.  Cross-section  data 

3.  Output  flux  energy  bin  limits 

4.  Weight  numbers  (for  each  re¬ 
gion) 

5.  Region-dependent  weights 


Problem  dependent:  (number,  E,  is  printed 
when  EZGEOM  is  run) 

Same  as  entry  1:  D  is  printed  (DATORG) 

J  (number  of  energy  bins)  +  1 
I  (number  of  regions) 

NRDW  (number  of  region-dependent 
weights) 
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Type  of  Data 


Number  of  Computer  Words  Used 


6.  Energy -dependent  weights 

7.  Composition  per  region 

8.  Absorption  per  region 

9.  Deaths  per  region 

10.  Degradations  per  region 

11.  Births  per  region 

12.  Absorption  per  element  per 
region 

13.  Cumulative  track-length  in  each 
region  and  energy  group 

14.  Cumulative  track-lengtiis 

15.  Cumulative  square  track-lengths 


(IDWT-NRDW)J  (IDWT=  total  number  of 
distinct  weights) 

I 

I 

I 

I 

I 

5«I 

J*I 

J*I 

J*I 


The  use  of  variable  dimensioning  allows  for  an  interplay  of  energy  bins,  regions, 
cross-section  data,  and  total  geometry  data.  As  long  as  the  total  number  of  words 
does  not  exceed  15,000,  one  may  run  an  ADONIS  problem. 

D  +  E  +  111  +  3JXI  +  J  +  1  +  NRDW  +  J(IDWT  -  NFJ)W)  ^  15,000 
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Department  of  the  Navy 

Washington,  D.C.  20360  1 

Commander 

U.  S.  Naval  Weapons  Laboratory 

Dahlgren,  Virginia  22148  2 


Commanding  Officer  and  Director 
U.  S.  Naval  Radiological  Defense 
Laboratory 
Attn:  W.  E.  Kreger 


San  Francisco  24,  California .  3 

AFWL  (SWOI) 

Kirtland  Air  Force  Base 

New  Mexico  87117  .  2 

AFCRL  (Dr.  Marcus  O’Day) 

L.  G.  Hans  com  Field 

Bedford,  Massachusetts  01731 .  1 

AFML  (MAAM) 

Wright-Patterson  Air  Force  Base 

Ohio  45433  2 

Director,  Project  RAND 
Department  of  the  Air  Force 
1700  Main  Street 

Santa  Monica,  California  90406  .  1 

National  Academy  of  Sciences 
2101  Constitution  Avenue,  N.W. 

Attn:  Mr.  R.  Park 

Washington  25,  D.C .  1 


Director 

National  Bureau  of  Standards 
Attn:  M.  Berger 

C.  Eisenhauer 
L.  V.  Spencer 
232  Dynamometer  Building 

Washington,  D.C.  20234  .  3 
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U.  S.  Public  Health  Service 
Washington,  D.C . 

U.  S.  Atomic  Energy  Commission 
Attn:  Document  Library 
Research  Division 

Washington,  D.C.  20545  . 

U.  S.  Atomic  Energy  Commission 
Director,  Civil  Effects  Test 
Group 

Attn:  Mr.  L.  J.  Deal 

Washington,  D.C.  20545  . 

U.  S.  Atomic  Energy  Commission 
Attn:  Dr.  Thomas  H.  Johnson 
Director 

Division  of  Research 
Washington  7,  D.C . 

U.  S.  Atomic  Energy  Commission 
Army  Reactors 

Division  of  Reactor  Development 
Washington,  D.C.  20545  . 

U.  S.  Atomic  Energy  Commission 
Technical  Information  Service 
Attn:  Dr.  H.  M.  Roth 
P.  O.  Box  62 

Oak  Ridge,  Tennessee . . 

University  of  California 
Lawrence  Radiation  Laboratory 
Attn:  Restricted  Data  Control  Officer 
P.  O.  Box  808 

Livermore,  California  94550  ... 

U.  S.  Atomic  Energy  Commission 
Los  Alamos  Scientific  Lab. 

Attn:  Technical  Library 
P.  O.  Box  1663 

Los  Alamos,  New  Mexico  87544  .  . 

Oak  Ridge  National  Laboratory 
Attn:  E.  P.  Blizzard 
Dr.  W.  Zobel 
F.  H.  Clark 
Dr.  J.  Auzier 

Oak  Ridge,  Tennessee . 


Argonne  National  Laboratory 
9700  South  Cass  Avenue 
Lemont,  Illinois  . 

Brookhaven  National  Laboratory 

Upton,  Long  Island 

New  York  . 

Aerojet -General  Corporation 
Azusa,  California  . 

Argonne  Cancer  Research  Hospital 
Attn:  Frances  J.  Skozen 

Editorial  Department 
950  E.  59th  Street 
Chicago  37,  Illinois  . 

Battelle  Memorial  Institute 
505  King  Avenue 
Attn:  Dr.  H.  W.  Russell 
Columbus  1,  Ohio  . 

General  Electric  Company 
Knolls  Atomic  Power  Laboratory 
Schenectady,  New  York . 

Lockheed  Aircraft  Corporation 
Lockheed  Nuclear  Products 
Attn:  M.  Burrell 
Georgia  Division 

Marietta,  Georgia  . 

Lockheed  Missiles  and  Space 
Division 

Technical  Information  Center 

3251  Hanover  Street 

Palo  Alto,  California . 

Space  Technology  Laboratories 
Attn:  Dr.  3.  Sussholz 
P.  O.  Bo.x  45565,  Airport  Station 
Los  Angeles  45,  California  .  .  . 

Technical  Operations,  Inc. 

Attn:  Dr.  Clark 
S.  Holland 
South  Avenue 

Burlington.  Massachusetts  .  .  . 


Westinghouse  Electric  Corporation 
Atomic  Power  Division 
Pittsburgh  30,  Pennsylvania 

Ottawa  University 
Department  of  Nuclear  Engineering 
Attn:  Dr.  L.  V.  Spencer 
Ottawa,  Kansas  . 

Australian  Group 
c/o  Military  Attache 
Australian  Embassy 
2001  Connecticut  Avenue,  N.W. 
Washington,  D.C.  20008 

The  Scientific  Information  Office 
Defense  Research  Staff 
British  Embassy 

31000  Massachusetts  Avenue,  N.W. 
Washington,  D.C.  20008 

Defense  Research  Member 
Canadian  Joint  Staff 
2450  Massachusetts  Avenue,  N.W. 
V/ashiJ.igton,  D.C.  20008 

Commanding  Officer 
Ballistic  Research  Laboratories 
Aberdeen  Proving  Ground 
Attn:  Mr.  Alfred  J.  Budka 
Aberdeen  Proving  Ground,  Md 


